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Abstract 

. We consider a Markov process X which is the solution of a stochastic differential equation driven by a Levy 

£Nj ■ process Z and an independent Wiener process W . Under some regularity conditions, including non-degeneracy 

of the diffusive and jump components of the process as well as smoothness of the Levy density of Z outside any 
neighborhood of the origin, we obtain a small-time second-order polynomial expansion for the tail distribution 
and the transition density of the process X. Our method of proof combines a recent regularizing technique 
for deriving the analog small-time expansions for a Levy process with some new tail and density estimates for 
^ C*j jump-diffusion processes with small jumps based on the theory of Malliavin calculus, flow of diffeomorphisms 



for SDEs, and time-reversibility. As an application, the leading term for out-of-the-money option prices in short 
maturity under a local jump-diffusion model is also derived. 



en ■ 1 Introduction 

> 

The small-time asymptotic behavior of the transition densities of Markov processes {X t (x)}t>o with deterministic 
. initial condition Xq(x) = x has been studied for a long time, with a certain focus to consider either purely- 
£f) ' continuous or purely-discontinuous processes. Starting from the problem of existence, there are several sets of 
sufficient conditions for the existence of the transition density of X t {x), denoted hereafter p t (-;x). A stream in this 
" direction is based on the machinery of Malliavin calculus, originally developed for continuous diffusions (see the 
monograph [35]) and, then, extended to Markov process with jumps (see the monograph [7]). This approach can 
also yield estimates of the transition density p t (-;x) in small time t. For purely-jump Markov processes, the key 
assumption is that the Levy measure of the process admits a smooth Levy density. The pioneer of this approach 
was Leandre [19], who obtained the first-order small-time asymptotic behavior of the transition density for fully 
supported Levy densities. This result was extended in |17j to the case where the point y cannot be reached with 



only one jump from x but rather with finitely many jumps, while |27| developed a method that can also be applied 
to Levy measures with a nonzero singular component (see also [35] and [T5J for other related results). 

The main result in []j|| states that, for y / 0, 

t 1 ^ ~i Pt ( x + y; ^ = 9 ^ X] y )' 

where g(x;y) is the so-called Levy density of the process X to be defined below (see (JSJ). Leandre's approach 
consisted of first separating the small jumps (say, those with sizes smaller than an e > 0) and the large jumps of 
the underlying Levy process, and then conditioning on the number of large jumps by time t. Malliavin's calculus 
was then applied to control the resulting density given that there is no large jump. For e > small enough, the 
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term when there is only one large jump was proved to be equivalent, up to a remainder of order o(t), to the term 
resulting from a model in which there is no small-jump component at all. Finally, the terms when there is more 
than one large jump were shown to be of order 0(t 2 ). 

Higher order expansions of the transition density of Markov processes with jumps have been considered quite 
recently and only for processes with finite jump activity (see, e.g., |35j ) or for Levy processes with possibly infinite 
jump-activity. We focus on the literature of the latter case due to its close connection to the present work. [32] 
was the first work to consider higher order expansions for the transition densities of Levy processes using Leandre's 
approach. Concretely, the following expansion for the transition densities {pt{y)}t>o of a Levy process {Z t }t>o was 
proposed therein: 

p t (y):=±F(Z t <y)= y^a n {y)-. +0(t N ), (y + 0, JV e N). (1) 
dy ^ n\ 

As in |19j . the idea was to justify that each higher order term (say, the term corresponding to k large jumps) can 
be replaced, up to a remainder of order 0(t N ), by the resulting density as if there were no small-jump component. 
However, this approach is able to produce the correct expressions for the higher order coefficients 02(2/), . . . only in 
the compound Poisson case (cf. [H]). The problem was subsequently resolved in [T3] (see Section 6 therein as well 
as |12j for a preliminary related result), using a new approach, under the assumption that the Levy density of the 
Levy process {Z t }t>o is sufficiently smooth and bounded outside any neighborhood of the origin. There are two key 
ideas in |12[ 113) . Firstly, instead of working directly with the transition densities, the following analog expansions 
for the tail probabilities were first obtained: 

nZt > y) = An ^n\ + (V>^ N ^ N), (2) 

n=l 

where sup 0<t<to \TZ t (y)\ < 00, for some to > 0. Secondly, by considering a smooth thresholding of the large jumps 
(so that the density of large jumps is smooth) and conditioning on the size of the first jump, it was possible to 
regularize the discontinuous functional l{z t >x} and, subsequently, proceed to use an iterated Dynkin's formula (see 
Section ^. 2l below for more information) to expand the resulting smooth moment functions E(/(Z t )) as a power series 
in t. ([T]) was then obtained by differentiation of ([2]), after justifying that the functions A n (y) and the remainder 
1Zt{y) were differentiable in y. 

The results and techniques described in the previous paragraph open the door to the study of higher order 
expansions for the transition densities of more general Markov models with infinite jump-activity. We take the 
analysis one step further and consider a jump-diffusion model with non-degenerate diffusion and jump components. 
Our analysis can also be applied to purely-discontinuous processes as in |19) . but we prefer to consider a "mixture 
model" due to its relevance in financial applications where compelling empirical evidence supports models containing 
both continuous and jump components (see Section|6]bclow for detailed references in this direction). More concretely, 
we consider the following Stochastic Differential Equations (SDE) driven by a Wiener process {Wt}t>o and an 
independent purcly-jump Levy process {Z t } t >o of bounded variation: 

X t (x) = x + j b(X u (x))du+ j a(X u (x))dW u + ]T j (X u -(x), AZ U ) . (3) 

uG(0,t]: AZ u7 tO 

Here, AZ U := Z u — Z u - := Z u — lim s /i t Z s denotes the jump Z at time u, and b, <r : K. — > M., 7:RxM- > M. are 
some suitable deterministic functions. The boundedness of the variation of Z is not essential for our analysis and 
is only imposed for the easiness of notation. 

As it will be evident from our work, an important difficulty to deal with the model (j3j) arises from the more 
complex interplay of the jump and continuous components. In particular, conditioning on the first "big jump" 
of {X s (x)} s <t leads us to consider the short-time expansions of the tail probability of a SDE with random initial 
value J, which creates important, albeit interesting, subtleties. More concretely, in the case of a Levy process (i.e. 
when b, a, and 7 above are state- independent), conditioning on the first big jump naturally leads to analyzing the 
small-time expansion of the tail probability P(A t e (x) + J > x + y), where {Xg(x)} stands for the "small jump" 
component of {X s (x)} (see the end of Section [2] for the terminology). This task is relatively simple to handle 
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since the smooth density of J "regularizes" the problem. In contrast, in the general local jump-diffusion model, 
conditioning on the first big jump leads to consider P(A t e (x + J) > x + y), a problem that does not allow a direct 
application of Dynkin's formula. Instead, to obtain the second-order expansion of the latter tail probability, we 
need to rely on smooth approximations of the tail probability building on the theoretical machinery of the flow of 
diffcomorphisms for SDEs and time-reversibility. 

Under certain regularity conditions on b, a and 7, as well as the Levy measure v of Z, we show the following 
second order expansion (as t — > 0) for the tail distribution of {X t (x)} t >o- 

t 2 

¥(X t (x)>x + y)=tA 1 (x;y) + -A 2 (x;y)+0(t 3 ), for x e R, y > 0. (4) 

The assumptions required for Q include boundedness and sufficient smoothness of the SDE's coefficients as well as 
non-degeneracy conditions on \d^j(x, £)| and |1 + d x "/(x, As in |19j . the key assumption on the Levy measure v 
of Z is that this admits a density h : K\{0} — > M. + that is bounded and sufficiently smooth outside any neighborhood 
of the origin. In that case, the leading term A\(x; y) depends only on the jump component of the process as follows 

A 1 (x; y) = v ({£ : j(x, Q > y}) = f h(Q<%. 

J{C---y(.x>0>v} 

The second order term ^2(2;; y) admits a more complex (but explicit) representation, which enables us, for instance, 
to precisely characterize the effects of the drift b and the diffusion a of the process in the likelihood of a "large" 
positive move (say, a move of size more than y) during a short time period t (see Remark 14.31 below for further 
details). 

Once the asymptotic expansion for tail distribution is obtained, we proceed to obtain a second order expansion 
for the transition density function pt(y; x). As expected from taking formal differentiation of the tail expansion (Q| 
with respect to y, the leading term of pt(x + y; x) is of the form tg{x;y) for y > 0, where g(x;y) is the so-called 
Levy density of the process {X t (x)} t > defined by 

9(x; y) := -|-f ({C : 7(*, > v}) (y > 0), (5) 
dy 

while the second order term takes the form — dyA^x; y) t 2 /2. One of the main subtleties here arises from attempting 
to control the density of X t (x) given that there is no "large" jump. To this end, we generalize the result in [TO] to 
the case where there is a non-degenerate diffusion component. Again, Malliavin calculus is proved to be the key 
tool for this task. 

Let us briefly make some remarks about the practical relevance of our results. Short-time asymptotics for the 
transition densities and distributions of Markov processes are important tools in many applications such as non- 
parametric estimation methods of the model under high-frequency sampling data and numerical approximations 
of functionals of the form $t(x) := K(4>(Xt(x)). In many of these applications, a certain discretization of the 
continuous-time object under study is needed and, in that case, short-time asymptotics are important not only in 
developing such discrete-time approximations but also to determine the rate of convergence of the discrete-time 
proxies to their continuous-time counterparts. 

As an instance of the applications referred to in the previous paragraph, a problem that has received a great deal 
of attention in the last few years is the study of small-time asymptotics for option prices and implied volatilities (see, 
e.g., PI], [TO], [U, 0. 03]> El, [31], H3, HI, [IS])- As a byproduct of the asymptotics for the tail distributions 
(HI, we derive here the leading term of the small-time expansion for the arbitrage- free prices of out-of-the-money 
European call options. Specifically, let {St}t>o be the stock price process and denote X t = log.S't for each t > 0. 
We assume that P is the option pricing measure and that under this measure the process {X t } t >o is of the form in 
(j3]). Then, we prove that 

Jim -E(S t - K)+ = / [S e^ - k) h(C)d(, (6) 

which extends the analog result for exponential Levy model (cf. [31] and [34]). A related paper is [21], where ([6]) 
was obtained for a wide class of multi-factor Levy Markov models under certain technical conditions (see Theorem 
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2.1 therein), including the requirement that lim t _s.o E(5t — K) + /t exists in the "out-of-the-money region" and some 
stringent integrability conditions on the Levy density h. 

The paper is organized as follows. In Sectional we introduced the model and the assumptions needed for our 
results. The probabilistic tools, such as the iterated Dynkin's formula as well as tail estimates for semimartingalcs 
with bounded jumps, are presented in Section [3l The main results of the paper are then stated in Sections [4] and [5l 
where the second order expansion for the tail distributions and the transition densities are obtained, respectively. 
The application of the expansion for the tail distribution to option pricing in local jump-diffusion financial models 
is presented in Section |6] The proofs of our main results as well as some preliminaries of Malliavin calculus on 
Wiener-Poisson spaces are given in several appendices. 



2 Setup, assumptions, and notation 

Throughout, C^ 1 (resp., C^°) represents the class of continuous (resp., bounded) functions with bounded and 
continuous partial derivatives of arbitrary order n > 1. We let {Zt}t>o be a driftless pure-jump Levy process of 
bounded variation with Levy measure v and {Wt}t>o be a Wiener process independent of Z, both of which are 
defined on a complete probability space (fi, T, P), equipped with the natural filtration (J-t)t>o generated by W and 
Z and augmented by all the null sets in T so that it satisfies the usual conditions (see, e.g., Chapter I in [30] ). 

As stated in the introduction, in this paper we consider the following local jump-diffusion model 

X t (x)=x + f b(X u (x))du+ f a(X u {x))dW u + V 7 (X u - (x), AZ U ) , (7) 
J ° J ° 0<u<t 

where AZ U := Z u — Z u - := Z u — \\va S/ At Z s denotes the jump Z at time u, and b, a : R — > R and 7 : R x R — > R are 
deterministic functions satisfying suitable conditions under which ([7} admits a unique solution. Typical sufficient 
conditions for ([7]) to be well-posed include linear growth and Lipschitz continuity of the coefficients b, a, and 7 (see, 
e.g., [31 Theorem 6.2.3], [251 Theorem 1.19]). It is convenient to write Z in terms of its Levy-Ito representation: 



Z t = V AZ U = f f (M(du,d0- 



0<u<t J ° jRo 

Here, Ro := R\{0} and M(du, d() := #{u > : (u, AZ U ) 6 dux d(} is the jump measure of the process Z, which is 
a Poisson random measure on R + x Ro with mean measure KM (du, d() = duv(dQ. We make the following standing 
assumptions about Z: 

(CI) The Levy measure v of Z has a C°°(R\{0}) strictly positive density h such that, for every e > 
and n > 0, 

(i) /(l A ICDMCR < 00, and (ii) sup |/i (n) (C)| < 00. (8) 

J |C|>B 

Remark 2.1 Condition fH-ij is equivalent to the condition that Z is of bounded variation. This condition is not 
essential for our results and is imposed only for the easiness of notation. We could have taken a general Levy process 
Z and replace the last term in ^ with a finite- variation pure-jump process plus a compensated Poisson integral as 
follows: 

X t (x)=x+ [ b(X u {x))du+ ( a(X u (x))dW u 
Jo Jo 

y(X u -(x),QM(du,d{)+ I I j(X u -(x),0M(du,d(), (9) 
C I > 1 Jo J\c\<i 

where M{du,dC ) ) is the compensated Poisson measure M{du,dC ) ) — v{dC ) )du. Condition is actually necessary 

for the tail probabilities of {X t (x)}t>a to admit an expansion in integer powers of time. Indeed, even in the simplest 
pure Levy case (X t (x) = Z t + x), it is possible to build examples where F(Z t > y) converges to at a fractional 
power oft in the absence of (see ]22p. 
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Throughout the paper, the generator 7 is assumed to satisfy the following conditions: 

(C2-a) •) e CpCR x R) and 7(2;, 0) = for all iel; 

(C2-b) There exists a constant 8 > such that \d^j(x, C)| > 5, for all 

Both of the previous conditions were also imposed in [E5]. Note that (C2-a) implies that, for any e > 0, there 
exists C e < 00 such that 

^7(^,0 



sup 



dx % 



< c £ \c\, (10) 



for all |C| < e and i > 0. Condition (C2-b) is imposed so that, for each the mapping £ — > 7(2;, £) admits an 

inverse function 7~ 1 (x, £) with bounded derivatives. Note that (C2-b) together with the continuity of d^(x, ()/d( 
implies that the mapping £ — > 7(2;, £) is either strictly increasing or decreasing for all x. 

We will also require the following boundedness and non-degeneracy conditions: 
(C3) The functions b(x) and v(x) := a 2 (x)/2 belong to C h °°(R). 
(C4) There exists a constant S > such that, for all i,(€R, 

<97(x, C) 



(i) 



Ox 



> 5, (ii) a{x) > 5. (11) 



Remark 2.2 Boundedness conditions of the type (C3) above are not restrictive in practice. Indeed, on one hand, 
extremely large values of b and a won't typically make sense in a particular financial or physical phenomenon 
in mind (e.g. a large volatility value a could hardly be justified financially). On the other hand, a stochastic 
model with arbitrary (but sufficiently regular) functions b and v could be closely approximated by a model with C£° 
functions b and v. The condition UlY i). which was also imposed in \ 19}/ . guarantees the a.s. existence of a flow 
<& St t(x) : R — > R, x — > X s j(x) of diffeomorphisms for all < s < t (cf. fl9f). where here {X s j(x)}t> s is defined as 
in ([7p but with initial condition X SiS (x) = x. 

As it is usually the case with Levy processes, we shall decompose Z into a compound Poisson process and 
a process with bounded jumps. More specifically, let (f> £ S C°°(R) be a truncation function such that l|f|> e < 
4>e(C) ^ l|ci>e/2 an d let {Zt(e)}t>o and {Z^.(e)}t>o be independent driftless Levy processes of bounded variation 
with respective Levy densities 

K(C) ■= MCMC) and Mc):=(i-&(C))M0- ( 12 ) 

Without loss of generality, we assume hereafter that 

Z t = Z t {e) + Z' t {e), for all t > 0. (13) 

The process Z'(e), that we referred to as the small-jump component of Z, is a pure-jump Levy process with jumps 
bounded by e. In contrast, the process Z(e), hereafter referred to as the big-jump component of Z, is a compound 
Poisson process with intensity of jumps A e := J (f) s (£)h(£)dC and jumps {Jf}%>i with probability density function 

MC>:=^>. (14) 

Throughout the paper, {A^ 6 } t >o denote the jump counting process of the compound Poisson process {Z t (e)}t>o 
and J := J e represent a random variable with density h s {Q. 

The following result, whose proof is presented in Appendix [D] will be needed in what follows. 

Lemma 2.3 Under the conditions (CI), (C2) and (C4), the following statements hold: 

1. Let 7(2, C) := 7(2,0 + z. Then, for each z £ R, the mapping f — > 7(2, Q (equiv. £ — > 7(2, £)) is invertible 
and its inverse 7 _1 (^,C) (resp. 7 _1 (z, Q) is C^ 1 (R x R). 
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2. Both j(z, J e ) and "f(z, J 6 ) admit densities in C™(KxR), denoted by I '(£; z) :~ T e (£; z) andT(£; z) := T e (£; z), 
respectively. Furthermore, they have the representation: 



r £ (c-,z) = h £ (T 1 (z,0) 



(*,r x M) 



r £ ((;z) = h £ ( 7 - 1 (z,0) 



(15) 



3. The mappings (z, () -> P (7(2, J £ ) > () and (z, () -> P (7(2, J e ) > C) are C b °°(R x R). 



^. TTie mapping z — > ti := z + 7(2, C) admits a density, denoted hereafter 7 (k, C)> belongs to Cf 1 (M. x M). 

Wc finish this section with the definition of two important processes. Throughout the paper, we let {X s (e, 0, x)} s>0 
be the solution of the SDE: 

X s (e 1 0),x):=x+ [ b(X u (e,Q,x))du+ f a (X u (e, 0, x)) dW u + ^ 7 (X u - (e, 0, x), A^( £ )) 

= x+ [ b(X u (e,V>,x))du+ [ o-(X u (e,(H,x))dWu+ [ [ J (X u -(e,$,x),() M'(du,d(), (16) 



where hereafter M' denotes the jump measure of the small-jump component Z'(e). The law of the process (fT6"|) can 
be interpreted as the law of {X s } < s <t conditioning on not having any "big" jumps during [0,i\. In other words, 
denoting the law of a process Y (resp. the conditional law of Y given an event B) by C(Y) (resp. C(Y\B)), we 
have that, for each fixed t > 0, 



C[{X s (x)} Q ^ t N? = 0)=£({X s (e,<b,x)} o < s ^ 
Similarly, for a collection of times < si < • • • < s„, let {X s (e, {si, . . . , s n }, a;)} s>0 be the solution of the SDE: 

X s (e, {si, . . . , s n },x) := x + / b(X u (e, {si, . . . , s„}, x))du + / a (X u (e, {si, . . . , s„}, x)) dW u 

Jo Jo 

+ l(X u -(e,{ Sl ,...,s n },x),AZ' u {e))+ ^ y(X ar (e,{si,...,s n },x),Jf), 



0<u<s 



i:Si<s 



where, as above, {Jf}i>i represents a random sample from the distribution </> E (C)^(C)^C/^e! independent of Z'. 
Then, denoting the arrival jumps of Z(e) by T\ < t 2 < . . . , it follows that 

C ( {X s (x)} 0<s<t N t =n,Ti = si,. . . ,r„ = s n ) = C ({X s (e, {si, . . .,s„},x)} 



fo<s<t 

The previous two processes will be needed in order to implement Leandre's approach in which the tail distribution 
¥(Xt(x) > x + y) is expanded in powers of time by conditioning on the number of jumps of Z(e) by time t. 

3 Probabilistic tools 

Throughout, C£(I) (resp. C^) denotes the class of functions having continuous and bounded derivatives of order 
< k < n in an open interval I C R (resp. in R). Also, ||g||oo = sup y \g(y)\- 



3.1 Uniform tail probability estimates 

The following general result will be important in the sequel. 

Proposition 3.1 Let M be a Poisson random measure on R + x Ro with mean measure MM(du, d£) = v{dC ) )dt and 
M be its compensated random measure. Let Y := Y^ be the solution of the SDE 



Y 



b(Y s )ds+ f a(Y s )dW s + f / 7(y s _,C)A/(ds,dC)• 
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Assume that b(x) and &(x) are uniformly bounded and */(x,Q is such that, for a constant S < 00, swp y |7(y,C)| < 
S(\£\ A 1), for v-a.e. C. In particular, the jumps of {Yt}t>o are bounded by S . Then there exists a constant C(S, k) 
depending on S and k, such that, for any fixed p > and all < t < 1, 



P<( sup 

k 0<s<* 



\Y S ( X} -x\ > 2pS^ < C(S,k)t p . 



Proof. Let V t be either J* a(Y s )dW s or J* J j(Y s _,z)M(ds,dz). It is clear that, in either case, V t is a mar- 
tingale with its jumps bounded by S. Moreover, its quadratic variation (V, V) is given by either J* a 2 (Y s )ds or 
Jo / 7 2 (^ s > O^idOds and, hence, in light of the boundedness of a and the 7, satisfies (V, V) t < kt, for some constant 
k. By equation (9) in [20], we have 



sup |Vs| > C } < 2exp 

0<s<t 



-AC + — kt(l + exp[A5] 



for all C, A > 0. 



(17) 



Now take C = 2pS and A = | \ogt\/2S. The rest of the proof is then clear. ■ 

As a direct corollary of the above proposition, we have the following tail probability estimate for the small-jump 
component {X t (e, 0, x)}t> of X defined in (TTJ 



Lemma 3.2 Fix any n > and a positive integer N. Then, under the Conditions (C2-C3) of Section^ there 
exist an e := e(N,n) > and C := C(N,rf) < 00 such that 



(1) For all t < 1, 



(2) For allt<l, 



sup 

0<e'<e,xGl 



\X t (e',Q,x)-x\ > rf) < Ct 



N 



(18) 



sup 



"{e 



X t (e' ,®,x)-x\ 



poo 

>s}ds= sup / P{|X f (e',0,x) -x\ > \ogs}ds < Ct N . 



Proof. The first statement is a special case of Proposition 13. 1[ which can be applied in light of the boundedness 
conditions (C3) as well as the condition (C2-a) that ensured that the components 



o-{X u {e'%,x))dW u: 



7 (X u -(e',Q),x),OM>{du,dO, 



are martingales with their quadratic variation bounded by kt, for some constant k > 0. Here, M' is the compensated 
Poisson measure of the small-jump component Z'(e'). To prove the second statement, we keep the notation of the 
proof of Proposition 13. II By (|17|) there exists a constant C > such that 



/•oo /• 00 

/ P{\X t (e,9,x) -x\ > log s}ds < C / exp 

Cr, 



A 2 

-A logs + yJW(l + exp[Ae]) 



ds 



(A - l)rr 



■ exp 



A 2 

—kt(l + exp[Ae] 



Now it suffices to take A = | log£|/2e and e = logr]/2N '. 



3.2 Iterated Dynkin's formula 

We now proceed to formally state a second-order iterated Dynkin's formula for the "small-time component" of X, 
{X t {e, 0, x)}t>o, defined in ([T5| . In general, the n-order iterated Dynkin's formula of a Markov process Y with 
initial condition Yq — x and generator Ly takes the generic form: 

n ~ 1 4-k MTL /*1 

= E J\ L ^ /(X) + Jn~lY. / (1 ~ ar ~ lE WW*)} da, (19) 
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where as usual L Y f = f and Lyf = Ly(Ly - /), n > 1. For the general jump-diffusion process ([9]), (fT9)) can be 
proved for n — 1 using Ito's formula (see |26[ Theorem 1.23]), while the general case can be proved by induction 
under somewhat strong conditions. In what follows, we only need the second-order formula and, thus, the subsequent 
result will suffice for our purposes. Below, L e denotes the infinitesimal generator for {X t (e, 0, x)} t >o which can be 
written as (c.f. [21)] ) 



L £ f(y)~VJ(y)+lJ(y), with 



(20) 



f"(y) + b(y)f'(y), Ief(y) ~ / (/(» + j(y, 0) - f{y)) MCK- 



Hereafter, h E (C) '■= (1 — ^(OMC) denotes the density of h truncated outside of a e-neighborhood of the origin. 
Let us remark that ([20]) is well-defined whenever / 6 C 6 2 in light of condition (C2-a). The proof of the following 
result is presented in Appendix [Pi 

Lemma 3.3 Under the Conditions (C1)-(C3) of Section^ the following assertions hold: 
1. For any function f in C\, 

Ef(X t {e,%,x)) = f{x)+t( EL £ f(X at (e,<D,x))da, (21) 

Jo 

where the remainder term is such that 



sup 

D<i<l,i£E 

s'<e 



EL e f{X at {e',%,x))da 



<K, 



(22) 



for a constant K < oo depending only on J \£\h e (Qd£, ||/^ fe '||oo> \\b^ ||oo, and \\v^ ||oo for k = 0, 1, 2. 
2. Additionally, if f G C*, then 

Ef (X t (e, 0, x)) = f(x) + tL s f(x) +t 2 f (1 - a)EL 2 J (X at (e, 0, x)) da, 

Jq 

where the remainder term is such that 



sup 

0<t<l,xGE 

e'<e 



(l-a)EL 2 e f(X at (e',H>,x))da 



<K, 



(23) 



(24) 



for a constant K < oo depending only on J |£|/i e (£)d£ , ||/^||ooj Wb^ Woo, and \\v^ k ' ||oo for k = 0, ... ,4. 



4 Second order expansion for the tail distributions 

We are ready to state our first main result; namely, we characterize the small-time behavior of the tail distribution 
of {X t (x)} t > : 

F t (x,y):=P(X t (x)>x + y), (y > 0). (25) 

As in [19] , the key idea is to use the decomposition (fT"3]) by conditioning on the number of "large" jumps occurring 
before time t. Concretely, recalling that {-/V t E } t > and A e := J E (C)MC)^C represent the jump counting process and 
the jump intensity of the large-jump component process {Z t (e)} t >o of Z, we have 

HMx) >x + y) = e-^ Y, nXt(x) > x + y\ JVf = n)±jj-. (26) 

n=0 
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The first term in (1261) (when n = 0) can be written as 

P(X t (x) > x + y\ Nl = 0) = F(X t (s, $,x)>x + y). 

In light of (fT8|) . this term can be made 0(t N ) for an arbitrarily large AT > 1, by taking e small enough. In order to 
deal with the other terms in (|26j) . we use the iterated Dynkin's formula introduced in Section [3.21 The following is 
the main result of this section (see Appendix [B] for the proof) . Below, h e and h £ denote the Levy densities defined 
in (fT2"|) . while g(x;y) denotes the so-called Levy density of the process {X t (x)} t >o defined by 

ff(x;y) :=-■£■ f h(Qd£ (27) 

oy J{(n(x,o>v} 

for y > 0. In light of Lemma 12.31 g admits the representation: 

g(x;y) = h{^~ x {x,y)) \ (d ( j) [x^~ l {x,y))\ 1 , 
where is the partial derivative of the function y(x, Q with respect to its second variable. 

Theorem 4.1 Let i£l and y > 0. Then, under the Conditions (C1-C4) of Section^ we have 

t 2 

F t (x,y) :=P{X t (x)>x + y)=tA 1 (x;y) + -A 2 (x;y) + 0(t 3 ), (28) 
as t — > 0, where A\(x]y) and A 2 (x;y) admit the following representations (for e > small enough): 



A^y) := / g{x;QdC= / h{QdQ, 
Jy J{~t(xX)>y} 

A 2 (x; y) := V(x; y) + Ji(x; y) + J 2 (x; y), 



with 



(29) 



g(x; ()d( + g(x; y) J + b{x + y)g{x\ y) 
a 2 (x) ( d 2 [°° d d \ 

+ —\d?J v !lin c)dC + 2 u;' ,Ln y) - 7J7' r: y) ) ' 

a(x + y) ( . . d . . .. . . A 

2 I + y)g^9[x; y) + 2<r (x + y)g{x; y) I 

// POO POO POO \ 

/ g(x;()dC+ _ g{x;QdC-2 g(x;C)d(\ h e (()d(, 

/• /*00 pOO p 

Mx;y)= g(x + j(xX);(')dChs(C)d(-2 g(x;()d( / h s (()d(. 

J Jy — t(x,0 Jy J 

Remark 4.2 Note that i/supp(z^) D {C : 7(2;, C) > 2/} = (s ^ * s 71 °^ possible to reach the level y from x with 
only one jump), then A\(x; y) = and P(X t (x) > x + y) = 0(t 2 ) as t — > 0. //, m addition, it is possible to reach 
the level y from x with two jumps, then J 2 {x\y) 7^ ; implying that ¥(X t (x) > x + y) decreases at the order oft 2 . 
These observations are consistent with the results in \ ITf and ]28jj . 

Remark 4.3 In the case that the generator j(x, £) does not depend on x, we get the following expansion for 
F(X t (x) >x + y): 

/>n ,> b(x)+b(x + y) , / a 2 (x) + a 2 (x + y) ... „ , ... . , A t 2 
9(()d( + V ' 2 V ^<?(y)* 2 - [ 2 9 + {X + V) {X + V)g{y) ) 2 

ff(CK- r 9(0dA h E (0d(t 2 



y-j(C) J y 

g(C)dCh E (()d( - 2 

y-7(() 
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The leading term in the above expression is determined by the jump component of the process and it has a natural 
interpretation: if within a very short time interval there is a "large" positive move (say, a move by more than 
y), this move must be due to a "large" jump. It is until the second term, when the diffusion and drift terms of 
the process X(x) appear. If, for instance, b and a are constants, the effect of a positive drift b > is to increase 
the probability of a "large" positive move of more than y by bg(y)t 2 (l + o(l)). Similarly, since typically g'(y) < 
when y > 0, the effect of a non-zero spot volatility a is to increase the probability of a "large" positive move by 
^| 5 '( y )|^(l + (l)). 



5 Expansion for the transition densities 

Our goal here is to obtain a second-order small-time approximation for the transition densities {pt(-~, x)}t>o of 
{X t (x)}t>Q. As it was done in the previous section, the idea is to work with the expansion (|26[) by first showing 
that each term there is differentiable with respect to y, and then determining their rates of convergence to as 
t — > 0. One of the main difficulties of this approach comes from controlling the term corresponding to no "large" 
jumps. As in the case of purely diffusion processes, Malliavin calculus is proved to be the key tool for this task. 
This analysis is presented in the following subsection before our main result is presented in Section [5.21 



5.1 Density estimates for SDE with bounded jumps 

In this part, we analyze the term corresponding to = 0: 

F{X t {x) > x + y\NI = O)=P(Jr t (e,0,aO >x + y). 

We will prove that, for any fixed positive integer N and r\ > 0, there exist an £o > and a constant C < oo (both 
only depending on N and rj) such that the density p t (-; e, 0, x) of X t (e, 0, x) satisfies 

sup p t (y;s,$,x) < Ct N , (30) 

\y-x\>V,E<co 

for all < t < 1. The estimate (|30| holds true for a larger class of SDE with bounded jumps. Concretely, throughout 
this section, we consider the more general equation 

Xf=x + f b(X*_)ds + f a(X*_)dW s + f [ 7 (X*_,0M(ds, dQ, (31) 
Jo Jo Jo J 

where, as before, M(ds, d() is a Poisson random measure on R + x K\{0} with mean measure p(ds, d() = dsx v{dC) 
and M = M — p is its compensated measure. As before, we assume that v admits a smooth density h with respect to 
the Lebesgue measure on R\{0}, but moreover, we now assume h is supported in a ball B(0, r) for some r € (0, oo). 
Note that the solution (X t (e, 0, x))t of (fT6)) perfectly fits within this framework. 

The following assumptions on the coefficient functions of (f3"Tj) are used in the sequel. Note that, when {Xf} is 
taken to be {X t (e, 0, x)}, these assumptions follow from the assumptions in Section [5] 



Assumption 5.1 The coefficients b(x),o~(x) and 7(2;, C) are 4 times differentiable. There exists a constant c > 
and a function k £ H2<p<oo ~^ p (y) such that: 



1. \d^b{x)\,\d^a{x)\, ^g?»7(*,0 



< c, for < n < 5. 



2. 



< c, for 1 < n + k < 5 and k > 1 . 



Malliavin calculus is the main tool to analyze the existence and smoothness of density for Xf . For the sake 
of completeness, we are presenting some necessary preliminaries of Malliavin calculus on Wiener-Poisson spaces 
in Appendix [3] following the approach in [7] . As described therein, there are different ways to define a Malliavin 
operator for such spaces. For our purposes, it suffices to consider the Malliavin operator corresponding to p = 
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(sec Section [A. 21 for the details). The intuitive explanation of p = is that when making perturbation of the sample 
path on the Wicner-Poisson space, we only perturb the Brownian path. 



Let us start by noting that Assumption 15.11 ensures that x —¥ Xf is a C -diffcomorphism with a continuous 
density (see, e.g., [7]). Furthermore, recalling that Hoo is the extended domain of the Malliavin operator, one has 
Xf £ Hoo and 

U t := T(X?,Xf) = (J a 2 (X*)J t (x)- 2 ds} J t (x) 2 . (32) 
Throughout this Section, we use the standard notation in which 

«,) - ±£. (33) 



Remark 5.2 Under the Condition (C4) of Section^ Jt{x) admits an inverse Y t := Jt{x) 1 , almost surely. Indeed, 
one can show that (cf. jTjj) 

dJ t (x) = 1 + d x b(Xf_)J t -(x)dt + d x o-{X*_)J t -{x)dW t 
+d x7 (X?_,0Jt-(x)M(dt,d0, 

while Y t = Jt(x)^ 1 satisfies an equation of the form: 

dY t = l + Yt-Dtdt + Y t -E t dW t + Y t _F t M(dt, dQ. 

Here Dt,Et and Ft are determined by b(x),o~(x),j(x,C) and Xf . As a consequence, together with our assumption 
on b, a and 7, one has 

E sup Jt(x) p , and E sup Jt{x)~ p < 00 
o<t<i o<t<i 

for all p > 1 . 

The main result of this section is Theorem 15 . 71 below. For this purpose, we state some preliminary known results. 
Let us start with the following integration by parts formula (the main ingredient for existence and smoothness of 
the density of Xf ), which is a special case of Lemma 4-14 in [7] together with the discussion of Chapter IV therein. 

Proposition 5.3 (Integration by Parts) For any f € C%°(R), there exists a random variable G t £ L p for all 

p£N, such that 

Ed x f(X?)=EG t U t - 2 f(X?). 

The following existence and regularity result for the density of a finite measure is well known (see, e.g., Theorem 
5.3 in [33]). 

Proposition 5.4 Let m be a finite measure supported in an open set O C K". Take any p > n. Suppose that there 
exists g = (gi, ...,g n ) G L p (m) such that 

f d J fdm= f fg s dm, /eC c °°(0), j = l,..,n. 

Then m has a bounded density function q £ Cb(O) satisfying 

IMloo < C\\g\\l P{m) m{Of- n l p . 

Here the constant C depends on n and p. 

The following lemma is where we use assumption (C4). It is the first ingredient in proving Theorem 15.71 
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Lemma 5.5 Recall Ut = T(X t ,X t ). Under the Condition (C4) of Section^ we have 

EU^ < Ct~P, 

for all p > 1 . 

Proof. The proof is a direct consequence of assumption (C4) and Remark 15.21 More precisely, 



U7 P = E- < — E 



(f*J s (x)-MXS) 2 ds) P " tP $ 2p ™io<s<tJs(x)-^ 



The proof is completed. ■ 

The last ingredient toward our main result of this section is the proposition below, which is just a restatement 
of Proposition 13.11 

Proposition 5.6 Let X t be the solution to equation \31\) . Recall that since the Levy density h is compactly supported 
and 7(x, z) is uniformly bounded, the jumps of X t are bounded by some constant S. There exists a constant C(S) 
depending on S , such that, for any fixed p > and all <t < 1, 

sup \Xf -x\> 2 P s\ < C(S)t p . 

0<s<t J 



Finally, we can state and prove our main result of this section. 



Theorem 5.7 Assume the Condition (C3) of Section^ is satisfied. Let {Xf} t >o be the solution to equation i31\) 
and denote by p t (y;x) the density of X£ . Fix r\ > and N > 0. Then, there exists r(rj,N) > such that if v is 
supported in £>(0, r) with r < r(j), N), we have for all < t < 1 

sup Vt {y;x)<C{r,,N)t N . 

\x—y\>r] 



Proof. For a fix t > 0, define a finite measure on R by 

m%(A)=F({X? eAr\B c (x,n)}), Aci, 

where B c (x,r) denotes the complement of closure of B(x,r). Thus, to prove our result it suffices to prove that 
ral admits a density that has the desired bound. To this end, for any smooth function / compactly supported in 
B c (x,rf), we have: 

/ (d x f)(y)mUdy) = m x f(xn - EG t U t - 2 f(xn 
Jr 

= J^E[G t U t - 2 \X? =y]mm?(dy), 

where the second equality follows from integration by parts. The rest of the proof follows from Proposition 15. 4[ 
Lemma 15.51 and Proposition 15.61 ■ 



5.2 Expansion for the transition density 

We are ready to state our main result of this section, namely, the second order expansion for the transition densities 
{pt(-', x)}t>o of the process {X t (x)}t>o m terms of the Levy density g(x; y) defined in (|27|) . The proof of the following 
result is presented in Appendix [C] 
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Theorem 5.8 Let x £ M and y > 0. Then, under the conditions and notation of Theorem ^. 1\ we have 




(34) 



as t — > 0, where a\(x; y) and 02(0:; y) admit the following representations (for e > small enough): 



ai(x;y) ■= g(x;y), 



a2(x;y) ■■= 3(x;y) + 9fi(x;y) + 3 2 (a;;y), 



3(a;;y) 



^(x;y) 



(35) 



9i(x;y) 



3 2 (x;y) 




(9 (x; y - j(x, 0) + g (x; 7(x + y, C) - x) <9„7(x + y, C) - 25(2; y)) h E (()d(, 
g(x + 1 (x,0;y-l(x,0)h s (Qd<:-2g(x;y) [ MCK- 



6 The first order term of the option price expansion 

In this section, we use our previous results to derive the leading term of the small-time expansion for option prices 
of out-of-the-money (OTM) European call options. This can be achieved by either the asymptotics of the tail 
distributions or the transition density. Given that the former requires less stringent conditions on the coefficients 
of the SDE, we choose the former approach. 

It is well-known by practitioners that the market implied volatility skewness is more pronounced as the expiration 
time approaches. Such a phenomenon indicates that a jump risk should be included into classical purely-continuous 
financial models (e.g. local volatility models and stochastic volatility models) to reproduce more accurately the 
implied volatility skews observed in short-term option prices. Moreover, further studies have shown that accurate 
modeling of the option market and asset prices requires a mixture of a continuous diffusive component and a jump 
component (see [T], [2], |J, [25], [S], and [23]). The study of small-time asymptotics of option prices and implied 
volatilities has grown significantly during the last decade, as it provides a convenient tool for testing various pricing 
models and calibrating parameters in each model (see, e.g., [16], [10], [14], [6], [11], [31], [34], [15], [24], [13]). In 
spite of the ample literature on the asymptotic behavior of the transition densities and option prices for either 
purely-continuous or purely-jump models, results on local jump-diffusion models are scarce. Our result in this 
section is thus a first attempt in this direction. 

Throughout this section, let {St}t>o be the stock price process and let X t = log 64 for each t > 0. We assume 
that P is the option pricing measure and that under this measure the process {X t }t>o is of the form in (|7]l. As usual, 
without loss of generality we assume that the risk-free interest rate r is 0. In particular, in order for St = expX t 
to be a Q- (local) martingale, we fix 



We assume that a and 7 are such that the Conditions (C1-C4) of Section [2] are satisfied. We also impose one extra 
condition (C5) for h{z) and 7(2;, z) in order to derive option price expansion, as we are working with exponential 
of jump-dissusion now. 

(C5) For each fixed x £ K there exists a constant M(x) > such that / e 3 ^^ x ' z ^h(z)dz < M(x) < 00. 

Note that this condition ensures immediately that b(x) above is well defined. 

By the Markov property of the system, it will suffice to compute a small-time expansion for 




Vi = E(S t - K)+ = E (e 




13 



In particular, using the well-known formula 

p oo 

EU1 {U>K} = KP{U > K}+ P{U > s}ds, 

Jk 

we can write 

poo poo 

E(e Xt - K)+ = / P{S t > s}ds = S P{X t -x> \ogs}ds, 
Jk Jk 

s 

where x = X = log So. Recall that 

°° (\ t\ n 

P(Xt -x>y) = e - A »* V V(X t - x > y\ N° = n)^-/-, (36) 

n=0 

where A £ := f cj) e (C)h(()d( is the jump intensity of {iV t e } t >o. We proceed as in Section|4j First, note that 

/■oo 
P{X t -x> \ogs}ds = Soe-^ih +I 2 + I 3 ), (37) 
K 

where 

/oo poo 
V{X t - x > \ogs\N~ = 0}ds = / ¥{X t (e, 0, x) - x > log s}ds, 

p OO 

I 2 = X e t / P{X t -x> logs|7V t e = l}ds, 

00 (\ f\n — 2 poo 

h = \lt 2 Y ^- Ll , / F{X t - x > log s\Nf = n}ds. 

n=2 J 

It is clear that I\ jt — >• as t — > by Lemma 13.21 We show that the same is true for I3 , which is the content of the 
following lemma. Its proof is given in Appendix |DJ 

Lemma 6.1 With the above notation, we have 



1 f°° 

sup — / P(\X t - a; I > log y\N^ = n)dy < 00. 
*,tero,xl n - Jo 



nGN,tG[0,l] m JO 

As a consequence, I^/t — > as t — > 0. 

Note that the above lemma actually implies that Ee'^'" 31 ' < 00 for all t G [0, 1). We are ready to state the main 
result of this section. 

Theorem 6.2 Let v t = E(S t — K) + be the price of a European call option with strike K . Under the Conditions 
(C1-C5), we have 



1 r / \ 

lim -v t = / (Soe***' - K) h(C)d(. 



Proof. We use the notation introduced in (|37p . Following a similar argument as in the proof of Lemma T6. 11 one 
can show that 

poo 

I sup P{X t - x > logs|iV t £ = l}ds < 00. (38) 

Also, it is clear that I\/t converges to when t approaches to by Lemma l3~2l Using the latter fact, (|38p . Lemma 
16. 1[ (|37l) , and dominated convergence theorem, we have 

G T poo POO 

lim — = lim = X e S lim / P{X t -x> \ags\Nf = l}ds = A,S / lim P{X t -x> \ogs\Nf = lids. 
t-s-0 t t-yO t t->o J K J k t->0 
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Next, using Theorem 14. 1[ it follows that 

pOC p OO p pOO 

\im— = S A x (x, log s)ds = S / / h(()d(ds = (5 e 7(2: - c) - K ) h(()d(, 

where we had used Fubini's theorem for the last equality above. ■ 

Remark 6.3 As a special case of our result, let 7(2;, £) = The model reduces to an exponential Levy model. The 
above first order asymptotics becomes to 

lim -^=/ (Soe< -K) + h(C)d(. 

t— >u 1 J-oo 

This recovers the well-known first-order asymptotic behavior for exponential Levy model (see, e.g., fSl]/ and \3^i ). 

A Preliminaries of Malliavin calculus on Wiener-Poisson spaces 

In this section, we follow [7] and introduce briefly the Malliavin calculus on Wicncr-Poisson space to our need. 
All the details can be found in [7]. We first introduce, in an abstract manner, the notion of Malliavin operator 
(infinitesimal generator of the Ornstein-Uhlenbeck semigroup). Fix a probability space (SI, J 7 , P). Let C 2 (M") be 
the space of all C 2 functions on K™ which, together with all their partial derivatives up to the second order, have 
at most polynomial growth. 

Definition A.l A Malliavin operator (L,1Z) is a linear operator L on a domain TZ C C\ p<00 ^ p , taking values in 
p| p<oo L p , and such that the following conditions hold true: 

1. If $ G 1Z n and F E C 2 (K"), then E TZ. 

2. T is the a-field generated by all functions in TZ. 

3. L is symmetric in L 2 , i.e. E($L*) = E(\E'L$) for all $,$eK. 
4- Define 

<J) = L($f) - - <JL$, for all $, * E TZ. 
Then T is nonnegative on TZ x TZ. 
5. For $ = ($\ $") S lZ n and F E C 2 (R"), one has 

" ap 1 " a2 p 

A.l Malliavin operator on Wiener-space 

We consider the Wiener space of continuous paths: 

H/ = (C([0,l],R),(J- t w/ )o<t<i,P w ) 

where: 

1. C([0, 1],K) is the space of continuous functions [0, 1] — > K; 

2. (W t ) t>0 is the coordinate process defined by W t (f) = / (t), f E C([0, 1],M); 

3. V w is the Wiener measure; 
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4. (J"^) <i<i is the (P^-completed) natural filtration of (W t ) 0<t<1 . 
Now, (fl w , T w ,V W ) is the canonical 1-dimensional Wiener space. Set 

5={<D = f(M/ (l ,..,W t J,FeC p 2 (nO<t 1 < 1}. 

For <E> given as above, define 

1 f)F 1 " f) 2 p 

lW * = - 2 E Wi Vr*>~; w t jw u + - £ ^(w tl ,..., w tn) (u a tj) . 

i—1 i ,j — 1 J 

It is well-known that (L w , 5 W ) is a Malliavin operator. Furthermore, for $ = F(W tl , Wi n ) and ^ = i?(W 5l , W 8m ) 

Tl Ti l O 771 O JT 

r ^*)=EEfc (m 1 ,...,^j^ T (^ 1 ,...,w 8ro )(^A Si ). 



»=1 3=1 



A. 2 Malliavin operator on Poisson space 

In this section, we consider the canonical probability space (H p , J rP ,V p ) of a prefixed Poisson random measure 
M(dt,d() on [0,1] x (R\{0}). Let p(dt,d() = dt x f(dC) be the intensity measure of M(dt,dQ and denote by 
M = M — fx its compensated Poisson random measure. Throughout our discussion, we also assume that v admits 
a smooth density h with respect to the Lebesgue measure on R, i.e. v{dQ = h(C)d(^. 

We denote by C^' 2 ([0, 1] x R\{0}) the set of all compactly supported, Borel functions / on [0, 1] x M\{0} that 
are of C 2 on R\{0} with /, d^f, d%f uniformly bounded. For such /, we write 



M(/) = J f(t,0M(dt,d(). 

Let p : R\{0} — > [0, oo) be an auxiliary function which is of class C\. Set 

S p = {$ = F(M(A), M(/fc)); F € C 2 (R fc ),/ i e C c b < 2 ([0, 1] x R\{0})}. 
Given any $ 6 <S P as above, define 

1 (9P / B h \ 

lP ® = 2 E d^( M (M> -' M (/*)) M f A/ t + 9,p5 c /, + p-£-a c /< J 

2—1 1 \ / 

In the above, stands for the Laplacian on M\{0}. 

Proposition A. 2 (L P ,S P ) defines a Malliavin operator. Moreover for 
$(M(/i), M(/ fe )) and * = H(M(hi), M(hj)), me We 

fiJP f)]-f 

r p ($, = E E fc( M W> — (M(/m), M^MipDcfiDchj). 

i—1 J — 1 * 

Proof. See Proposition 9-3 in [7J. ■ 

Remark A. 3 .As we can see from the above definition, the Malliavin operator for Poisson space is not uniquely 
defined. Indeed, each different choice of p gives a different operator. More remarks on this point are given at the 
end of this section below. 
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A. 3 Malliavin operator on Wiener-Poisson space 

Keep notations as in previous sections. We define a probability space (fi, J 7 , P) by 

(0,J",P) = {Vt w ,F W ,F W )® {Q P ,T P ,V P ). 

Clearly (O, T, P) is the canonical Wiener-Poisson space. Then define (L, S) to be the direct product of (L W ,S W ) 
and (L P ,S P ). More precisely, 

1. S = {$ = $„; *i, * m ); for all 6 e S p ,F e C 2 p }. 

2. For $ as given in (1), and w = (w 147 ,^) e f2, 

L$(a>) = i w/ $(-,w p ) + L p <5>{uj w , •). 

It can be shown that (L,S) is a Malliavin operator on (H,T,F). A few properties of (L,S) are listed below. 
Proposition A. 4 Let <f> e S w and >te5 p , we have 

1. L<P = L W <S>, and $) = T W (<S>, $). 

2. = L p #, and r(*, *) = r p (*, 

3. r($,*) = o. 

Finally, we extend the L to a larger domain that is suitable to work on. For this purpose, set 

\Mh p = + + ||r(*, *) 1/2 ||i^. 

and denote by H p the closure of <S under j| ■ \\h p - Also set 

Hoc = P) -Hp- 

2<p<oo 

Proposition A. 5 The operator {L,!!^) is a Malliavin operator. 
Proof. This is the content of Section 8-b in [7]. ■ 

Remark A. 6 As before, a different choice of the auxiliary function p in the definition for L p results in a different 
Malliavin operator L on the Wiener-Poisson space. Hence, the extended domain Hrx, may also be different. 

Remark A. 7 Consider the following stochastic differential equation with regular enough coefficients 

Xf=x + [ b(X*_)ds + [ a(X*_)dW s + f f j(Xf_,QM(ds,dO. 
Jo Jo Jo J 

To make use of the power of Malliavin calculus to analyze the above SDE, one needs to impose further conditions 
on p so that Xf is in the domain Hoc. Those further conditions are generally determined by the support of the 
intensity measure v and are needed to ensure that some approximation procedure can work out. For more details of 
the role that p plays, see Section 9-1 and Theorem 10-3 in It is important to note that, among various choices 
of p, one can always choose p = 0, as we impose the non- degeneracy Assumption (C3). The intuitive explanation 
of p = is that when making perturbation of the sample path on the Wiener-Poisson space, we only perturb the 
Brownian path. 



17 



B Proof of the tail distribution expansion 



The proof of Thcorem l4.1l is decomposed into three steps described in the following three subsections. For future use 
in obtaining the expansion for the transition densities, we will write explicitly the reminder terms when applying 
Dynkin's formula (|23[) or in any other type of approximation. 



B.l Key lemma to control the tail of the process with one large jump 

The following result will allow us to obtain the second-order expansion for the process with one large jump. 
Lemma B.l Under the setting and conditions (C1-C4) of Section^ 

P (X t (e, 0, z + 7 (z, J)) > 0) = H (z; 0) + tH^z; 0) + t 2 U t (z; 0), 

for any z, S M, where 

H {z;i9) :=P( 7 (z,J) + ;z>0), Fi(z;0) :=£>(*; 0) + /(*; 0), 
D{z; 0) := f (0; z)b{d) - d*f (0; z)v(4) - f (0; z)v'{&), 

I(z; 0) := / [P (z + j(z, J) > 7(0, 0) - P (* + l(z, J) > 0)] h e (()d{, 



(39) 



and, for e > small enough, 



lim sup sup 

t->0 zGR 



H\(z]'d) < oo, sup |-ffi(z; 0)| < oo. 



z.<9 



The idea to obtain Q39p consists of approximating the function l{x t (e.0, 2+7(2, J))>$} a smooth sequence of functions 
fs(X t (e, 0, z + j(z, J))), (5 > 0. Concretely, we let 



fs(w) 



ips(u)du, 



where <ps(x) := 6 1 f{8 1 x) for a density function <p £ C°° with supp^ C [—1, 1]. In particular, as 6 — > 0, 

f s (w) 1{ W >&}, and y g(w)f' s (w)dw = J g{w)tp s (w - tf)dw -> g(0), (40) 

whenever (7 is bounded and continuous at 0. It is then natural to apply Dynkin's formula to fs(X t (e, 0, z + ^(z, J))) 
and show that each of the resulting terms is convergent when 6 — > 0. The following result, whose proof is presented 
in Appendix Q21 is needed to formalize the last step. 

Lemma B.2 Let T(-;z) be the density of the random variable z + 7(2, J) and let pt (•;£, 0, C) be the density of 
Xt(e, 0,C)- Then, under the Conditions (C1-C4) of Section^ there exists an e > small enough such that for 
any compact set K C M, 



lim sup sup sup 

t-s-0 zSRijSif 



< OO, 



k > 0. 



(41) 



We are now in position to show (f3l)|). Proof of Lemma IB. 11 Throughout, <9 y7 and <9^ 7 will denote the 
partial derivatives of 7(2/, £) with respect to its first and second arguments, respectively. By dominated convergence 
theorem, we have 

(42) 



Note that 



' (X t {e, 0, z + 7(3, J)) > 0) = limE/ 5 (X t (e, 0, z + 7(2, J))). 

<5J,0 



E/ 5 (X t (e,0,z + 7 (z,J)))= / r(C;z)E/ a (X t (£,0,C)K 



(43) 
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and, thus, applying Dynkin's formula (|24[) to the expectation in the above integral gives us 

Ef s (X t (e, 0, z + 7 (z, J))) = J f(C; z)f s (()dC + 1 J f (C; z)L e f s (()dt (44) 

+ i 2 / f(C;z) / (1 - a)E(L e ) 2 f s (X at (e 7 tbX))dad(;. (45) 



J JO 

We analyze the limit of each of the three terms on the right-hand side of the previous equation. By dominated 
convergence theorem, the leading term of (|42[) is given by 

H Q (z;#) :=lim J ?((; z)f s (()d( = J f (£; z)I[# 00 - ) (()d( = P{j(z, J) + z > •&}. 

To compute the limit of the subsequent term, recall L e fs = T> e fg + I £ fs with V e and X e defined as in (|2T)|) . Then, 
the term of order t has two contribution: 



A S := / r(C;*)2? e /i(C)dC. B S := / r(C; z)X e f s (Od(. 
Using that fg(() — <fs(C ~ an d by integration by parts, it follows that 

A S = J (f(C,z)b(Q-d ( f(C;z)v(0-r(C;z)v'(0) MC-#)dC. 

Applying (gUJ* and Lemma 1231(2). 

lim A s = f (i?; z)6(i?) - d^T(i); z)v(i)) - f (#; z)v'(d). 

(54-0 

We now analyze the limit of the second term Bs. Applying (PT|) given in the Appendix [D] below, we can write Bs 

as 



Since 



Bs= I f' 5 (w) r(7j; z)d V h £ (()d(dw = <p s (w-<ff) / / T(r);z)drih £ (0d(dw. (46) 

f^; z)^MCR = - / [ r(y(wXP)\z)(da)( W> (WPhe(0(d( 

where we used ([80]) in the second equality, conditions (C2) and (C4) imply that the factor multiplying (ps(w — $) 
in (|46|) is bounded and continuous in w and, thus, 

\imB s = f [ f( V ;z)d7 1 h e (C)dC=:B (z;d). 

Next, recalling that T(C', z) is the density of J := z + j(z, J), 

B (z; 0) = /" (P (* + •/) > 7(*> 0) - P (« + 7(«, J) > <>)) MCR- 
Putting together the previous two limits, we obtain the term of order t: 

H^z; ??) := lim f f (C; z)L e f 5 (()dC = D(z; t?) + 7(z; 0), 
with D(.z; i?) and J(z; i?) given as in the statement of Lemma. 
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Finally, we estimate the remainder term 



1 

2 , 



Kt(z; 0) := lim / r(C; z) (1 - a)E(L e ) 2 f s (X at (e, 0, ())dad( (47) 

S i° J JO 

and show that this is uniformly bounded for t small enough. Let 1Zt(z;$;8,e) be the expression following lim^o 
and note that 

%t(z;#,6,e) = Jt{C;z)J o (1 - a)E(V £ ) 2 f s (X at (e, 0, ())dad<; 

+ J f(C;^ (1 - a)E(l £ ) 2 f s (X at (e,(bX))dad( (48) 
+ fr(C;z)J {l-a)KE s V £ f s {X at (e,®,Q)dad<; 
r(C; z)d( f (1 - a)EV e 2efs(X at (e, 0, Q)dad(. 



We shall use Lemmas IB. 21 and ID. 21 to deal with the four terms on the right-hand side of the previous equation. For 

—(2) 

simplicity, we only give the details for second term, that we denote hereafter J 4 (#;<$, e, z). The other terms can 
similarly be handled. Using (|W)) - (p?l"j) and Fubini, 

% 2) (z;#,S,e)= [ f(C;^ (!-«)/ %fs(w)p at (w; e, 0, ()dwdad( 
F((;z) / (1-a) / X s fs(w)X E p a t(w;e,<H,C)dwdadC 



r(C;z) / (1-a) / /» / / l £ p at (me,®,()dr)h e (OdCdwdad( 

JO J J Jy(w,Q 

pi p pw P 

f' s (w) (1-a) / T(C,z)l £Pat ( m eJX)d(dr,h e (C)d(dadw. 
Jo J Jj(w,C) J 

Using Fubini, it follows that 

?(C; z)X e p at (r); e, 0, Qd( = %Pat{m z , e), where p t (m z, e) : = / f (C; z)p t (ty e, 0, C) d(. 
Then, using the identity 

% 9 (y) = -//V(i( > ,Cffl(»<7)(,,CW ITI5 || M « 
-rt)/»0,fl 1 + M ^ M «, 



we have 



lf ) (z\$,8,e) = I f' s (w) / (1-a) / / T e p at (r); z,e)dr)h £ (()dCdadw 

Jo J J>y(w,() 

fs{w)I ( t 2,X) {w\z,£)dw+ / f' s (w)I^ 2 \w;z,e)dw, 
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where 



t(2,1)/ v f\, , f f W u , > ( {dyimjh Q,OMO ^ 7 r r.,r, ,. Q , 

J t 0;z,e):=- / (1-a) / / _p at (v,z,e) / " — -— — 7 r- 7 ^dC,dr]h e (QdC,da (49) 



7(to,C) 



ir\ W ;z,e):=- (1-a) / / / ^(^C^^.e))^)^^^ 

JO J Jj(w,() J JO 

-d(dr]h £ (C,)d(da. 



1 + (^7X7(77, 0,0 

We next show that each l[ 2,l \w, z,e) is uniformly bounded in w and z for 7; small enough. Indeed, using (fTUf and 

cih), 

|4 2,1) (w;z,e)| <C / (1-a) / / p at (mz,e)dr]h s (0dCda 

JO J J~/(wX)Aw 

<C [ (1-a) / [ p^jiw^^z^^da^w^m^rheiOdCda 
Jo J Jo 

<C f f f (l-a)p at (%w,(0)iz,e)dadl3(h £ {Od{, 
J Jo Jo 

for some constants C (that may change from line to line). Using (|4"Tj) with fc = 0, it follows that, for e,t > small 
enough, 

sup |I t (2,1) (w;;z,e)| < 00. (50) 

Similarly, 

|4 2 ' 2 W,e)l <C / (1-a) I / / |#, t (7fo, Ci8)5«,e)I^C^(C)dC*(C)dCda 

JO J Jj(w,QAw J Jo 

<C f (1-a) [ f [ /" U^CtCtC*", C^), C^)J ^. e) IdjSC ^ (C)<iC^ Cfc« (OdCdQf, 
Jo J Jo J Jo 

which, again using (|4"Tj) with k = 1, leads us to 

sup \I ( t 2 ' 2 \w;z,e)\ < 00, (51) 

z£R,?u£|supp/i 

for e, t > small enough. Due to (|50ti51p and the continuity of the functions ^(S, w; z, t) in 

lim/ t (2) (z; 0; 6, e) = lf' X) (§; z, e) + I t (2 ' 2) (0; z, e), (52) 
0^0 

which is uniformly bounded in z for any fixed 1? and < t < t Q with i > small enough. ■ 
B.2 The leading term 

In order to determine the leading term of (|25l) . we analyze the second term in (|26l) corresponding to n = 1 (only 
one "large" jump). By conditioning on the time of the jump (necessarily uniformly distributed on [0,t]), 

1 r* 

F(X t (x)>x + y\N? = 1) = - F(X t (e,{s},x)>x + y)ds. (53) 



* Jo 

Conditioning on T s - , 

F(X t (e, {s}, x) > x + y) = E (G t _ s (X s - (s, 0, x))) = E (G t _ s (X s (e, 0, x))) , (54) 
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where 



G t (z) :=G t {z;x,y) := P [X t (e, 0, z + 7 (z, J)) > x + y] 



(55) 



Using Lemma |B. 11 



¥{X t {e, {s}, x) > x + y) = EH (X s {e, 0, x); x + y) + (t- s)EH 1 (X s (e, 0, x);x + y) (56) 

+ (t - s) 2 E1ll_ s {X s {e,®,x);x,y), 

where TZl(w; x, y) := TZt(w; x + y). Next, we apply the Dynkin's formula (f2"3]l to EH {X s (e, 0, x); x + y), which is 
valid since H (z; x + y)=P (7(2, J) + z > x + y) is Cj* in light of Lemma CH(3). By 



lH Q (X s (e,Q),x);x + y) = H ,o(x;y) + sH ,i(x;y) + s TZ s {x;y) 



(57) 



where 



H 0fi (x; y) := H (x; x + y)=P [7(1, J) > y] , 

tt ( \ t tt 1 , \ u ,dH (z;x + y) 
Ho,i\x\y) := L £ H Q (x;x + y) = b{x)- 



a 2 (x) d 2 H (z;x + y) 



dz 



+ / (H (x + 1 (xX);x + y)- H (x;x + y))h £ {()d( 



dz 2 



b(x 



dP [ 7 (z, J) > x + y] 



dz 



a 2 (x) d 2 P [ 7 (z, J) > x + y] 



9z 2 



+ #o,i(a:; y), 



(58) 



Hj(x;y) := / (1 - a)E^ir (x as (e, 0, ar); x + y)cfa, 
with H 0i i(x;y) given by 

H ,i(x;y)= I (P [ 7 (x + 7(2;, (), J) >y- j(x, Q] - P [j{x, J) > y}) h £ (()d(. 



(59) 



Note that sup s<l x \TZ 2 (x; y)\ < 00 in light of ([23)1 and, also, by writing P [ 7 (z, J) > x + y] = P [7(2, J) > x + y — z] 
G(z, x + y — z) with G(x, y) = P ( 7 (x, J) > y), we get 



dP [j(z, J) > x + y] 



dz 

d 2 P [%z, J) > x + y] 
dz 2 



OP [ 7 (x, J) > y] 
dx 



+ T(y;x), 



S 2 P[ 7 (x,J) >y] , cT(y;x) cT(y;x) 



<9x 2 



+ 2 



<9x 



dy 



(60) 
(61) 



Plugging (|57p in (|56p and recalling from Lemma IB . II that the second and third terms on the right hand side of 
([56)) are bounded for t small enough, P(X t (e, {s}, x) > x + y) = P [j(x, J) > y] + 0(t), which can then be plugged 
in ([53]) to get 

P(X t (x) > x + y\ iV t £ = 1) = P [ 7 (x, J) > y] + 0(t). 

Finally, (|26]l can be written as 

P(X t (x) > x + y) = e-^H\ £ P [ 7 (x, J) > y] + 0(i 2 ) = t J l {l(xX) > y} h{QdC, + 0{t 2 ), 
where the second equality above is valid for e > small enough. 
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B.3 Second order term 



In addition to (f57|) . we also need to consider the leading terms in the term EHi(X s (e, 0, x); x + y) of ([56| and the 
term P(X t (x) >x + y\N£ = 2) of ([26]). Let us first show that z ->■ H^z; x + y) is C b 2 . To this end, let 

AC(C; i, y, ^) := F [z + 7 (z, J) > 7(x + y, Q] - P [« + 7(«, J) > a; + y], 

and recall that 

ffi(z;x + y) = ?(x + y;z)b(x + y) - (<%F)(x + y; z)v(x + y)- f(x + y;z)v'(x + y) + K,((;x,y,z)h e (C)d(, 



where d^T and d z T denote the partial derivatives of the density r(£;z). Obviously, the first three terms on the 
right-hand side of the previous expression are C% in light of Lemma I2~31 -f2). Hence, for the derivative d z Hi(z;x + y) 
to exist, it suffices to show that d z K,(Q] x, y, z) exists and that 

dK.(C;x,y,z) 



sup 



dz 



(62) 



for any |C| < £ and some constant C < oo. Wrtting /C(C; x, y, z) = fZ x+y q r(»7; z)dr\ and using that f (?j; z) e C£°, 
we have 

3K{t;x,y,z) f x+y dT( V ;z) 



dz 



dz 



-d?]. 



Therefore, 



\d z lC{Q;x,y,z)\ < sup d,T{r,;z) / |d c 7(x + y, C/3)|d/3|C| < sup d z T(r,;z) \dtf(x + y,ri)\ |C|, 

Tj.z Jo \x,y,z,T] J 

in light of Lemma 12.31 We can similarly prove that d z Hi(z; x,y) exists and is bounded. Using (|21[) and that 

f(C;z) =T(C-z;z), we get 

Effipf^e, 0,x); x, y) = H lfi (x, y) + s1l 3 s {x; y), 



(63) 



where H lfi (x\y) := H\{x;x + y) = T> lf0 (x;y) + Hi fi (x;y) with 

T>ifi{x;y) := T(y; x)b(x + y) - (d ( T)(y;x)v(x + y) - L(y; x)v'(x + y) 
H lfi {x; y) ■= I (P [x + j{x, J) > 7(x + y, 0] -P[x + y(x, J) > x + y]) MCK, 

K 3 s (x;y) 



[ WL e H 1 {X a8 {e,$,x);x + y)da = 0{l), as s -> 0. 
Jo 



In order to handle P(X t (x) > x + y| iVf = 2), we again condition on the times of the jumps, which are necessarily 
distributed as the order statistics of two independent uniform [0,i] random variables. Concretely, 



\X t (x)>x + y\NI =2) = ^ 



Js 



\X t {e, {si, s 2 }, x) > x + y)ds 2 dsi 



(64) 



Next, we determine the leading term of P(-Xt(e, {si, S2}, x) > x + y). By conditioning on T s - , 

P(X t (e,{ Sl ,s 2 },x) > x + y) =E(G t „ S2 (X S2 (e, {si},x))) , 

where, by Lemma IB. 11 

G t (z) = P [X t (e, 0, z + j(z, J)) > x + y] = H (z- x + y) + tflifo x + y) + t 2 TZ t (z; x + y). 



(65) 
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Then, for e > and t small enough, 

P(X t {e, {si, s 2 }, x) > x + y) = E (H {X S2 (e, { Sl }, x); x + y)) (66) 

+ (t - s 2 )EU^_ 82 (X S2 (e, {si}, x); x, y) , 
with Ht( z i x i v) := H\(z; x + y) + tTZt(z; x + y). 
Again, conditioning on J- - 

E (H a (X S2 (e,{si},x);x + y)) = E (X + yj), 

where 

G t {z- x + y) := EH Q (X t (e, 0, z + 7 (z, J)); a; + y) . 
Since z — > H (z; x + y) = P(z + -f(z, J) > x + y) is C£° by Lemma [231 - (3). one can apply (|2"Tj) to deduce 



G t (z;x + y)= / r(C; z)EiJ (X t (e, 0, C)i x + y) d( = I T((; z)H ({;x + y) d£ + HR*{z;x,y) 
=:H 2 (z;x + y)+tTZ e t (z;x,y), 
where, denoting two independent copies of J by J\, J 2 , 

H 2 (z; x + y):=P(z + j(z, Ji) + j(z + 7 (z, Ji), J 2 ) > x + y) , 

K 6 t (z;x,y):= Jr((;z)J^ EL £ H Q (X at (s, 0, (); x + y) dad(. 

Therefore, 

P(X t {e, {si, « a }, x) > £ + y) = E (7J 2 (X Sl (e, 0, x); x + y)) + (s 2 - s x )E^ 2 _ sl (X Sl (e, 0, x); x, y) 

+ (i - s 2 )Eftf_ S2 (X S2 (e,{si} ,x);x,y) . 
Applying again (|21l) to the first term on the right-hand side of the previous equation, 



P{X t (e, {si, s 2 }, x)>x + y) = H 2fi (x; y) + Sl TZ 5 Sl (x; y) (67) 

+ (s 2 - Sl )EK%_ Si (X S1 (e, 0, x); x, y) 
+ (t- s 2 )E7^_ S2 (X S2 {e,{ Sl } : x);x,y) , 



where 

H 2 , (x;y) := i? 2 (x;x + y) = P( 7 (a, Ji) + 7 (x + 7 (x, J x ), J 2 ) > y) 



/ EL £ H 2 (X 

< V .S' 1 

(e, 0,x);x + y)da 

Jo 



Therefore, we conclude that 

P (Xt(s) > x + y|A^ e = 2) = ff 2 , (x; y) + 0(t). (68) 

In light of (|5^| . (|56p - ([55|) , (jB"3")) . and we have the following second order decomposition of the tail distribution 

P(X t (x) >x + y): 

A t 2 

P{X t {x) > x + y) = e-^XetHoAx; y) + (Ho,i{ x > v) + #i,ofo v)) 

+ e- x ° t ^H 2 , (x;y) + O(t 3 ) 



= X e tH 0t o(x;y) + y|A e [H ,i(x;y) + Hi )0 (x;y)] 
+ \l [H 2 , a (x; y) - 2H Q , (x; y)\ } + 0(t 3 ). 
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The expressions in (|29|) follows from the fact that, for e small enough, 

poo p poo 

A £ P[ 7 (x,J) >y] = / \ e T s ((;x)d{= / MCR = / <?(*;CK- (69) 

J y J{C:7(z.C)>2/} 

Thus, for fixed i£l and y > and e small enough, 

\ s T e ((;x) = g(x;Q. (70) 
Using the two previous relationships, it is not hard to see that 



X s H ,o(x;y) = j g(x;C)dC; A e H ,i(x; y) + #1,0(2:; y) = Ji(x\y), 

Jy 

A £ [H ,i(x;y) +H h0 (x;y)] = V{x;y) + Ji(x;y), \ 2 e [H 2fi (x;y) - 2H , (x;y)} = J 2 (x;y). 
This concludes the result of Theorem 14.11 



C Proof of the expansion for the transition densities 

The following result will allow us to control the higher order terms (see Appendix iDl for a proof): 
Lemma C.l Let 

H t (x,y) := e"*-* V ¥(X t (x) >x + y\N? = „)M^. (71) 

Then, under the conditions of Theorem \5.8[ there exists e > small enough as well as to := to(s) > and 
B = B(e) < oo such that, for any < t < to, 

\d y Kt(x,y)\ < Bt 3 . 

Proof of Theorem 15.81 Let us consider the terms corresponding to one and two "large" jumps in (|26|) . From 
Ol. (1561). (IBTl). and (l63l). it follows that 



V(X t (x)>x + y\N^ = l) = H 0fi (x;y) + t[H 0tl (x;y)+H lfi (x;y)] (72) 
1 t* 

+ - / {s 2 Tl 2 s (x;y) + (t - s)sn 3 s (x;y) + (t - s) 2 E1ll_ s (X s (E,^xy,x,y)} ds. 
1 Jo 

Similarly, from (jS4"|). and (joTf . we have 

P(X t (x) > x + y\N e t = 2) = H 2 ,o{x; y) 

Sl K 5 Sl (x; y) + (s 2 - Sl )ETl 6 S2 _ Sl (X Sl (e, 0, x); x, y) (73) 

+ (t— s 2 )Ellt_ S2 (X S2 (e,{s 1 },x);x,y) \ds 2 dsi. 



2 

~~2 



f- 



Js 



Equations (|72j|73[) show that in order for the derivatives 

ai(x;y) := g-F(X t (x) >x + y\N? = 1), a 2 (x;y) := -P(l t ^) > x + y\N; = 2) 

to exist, it suffices that the partial derivatives of the functions Hij(x;y) with respect to y exist and also that the 
partial derivatives of the two types of functions, 1Z\{x;y) with i = 2,3,5 and lZ J t (w;x,y) with j = 1,4,6, with 
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respect to y exist and arc uniformly bounded onw£l and on a neighborhood of y. Furthermore, under the later 
boundedness property, we will then obtain that 



- , s dH 0fi (x;y) t 

«i(*;») = — — + 2 



dH 0A (x;y) dH h0 (x;y) 
dy dy 



0(t 2 ), (t^O), (74) 



a 2 (x;y)= dH2 f* ;y) +0(t), (t ^ 0). (75) 



Note that (|74ti75)) suffices to conclude (jMJ) in light of ([2l)]l. Theorem [S~71 and Lemma ICTTl 
(1) Differentiability of Hij(x;y): 

The differentiability of the functions Hij(x;y) essentially follows from Lemma \2. 3 1 Indeed, Lemma [2~3l -f2) implies 
that d y Ho,o(x;y) = —T(y;x) as well as 

a 2 (x) ( d 2 T{y;x) n d 2 T(y;x) d 2 T(y;x) \ ( dT(y;x) , dT(y;x) 

d y H s(x; y) := — — — — + 2 —5— + b(x) 



dx 2 dydx dy 2 J \ dx dy 

+ I (T(y;x)-T(y- 1 (x,0;x + j(x,0)h s (OdC- 



h 2 A^ y) = s- / p (7(s + 7(*, CO, Ji)>y- i{x, a)) MC1R1 



Similarly, we have that 

d y H lfi (x;y) := d y V lfl {a 
To compute d y H2fi(x;y), note that 

■^-H 2 ,o{x;y / 
dy dy 

|- r r(C 2 ;x + 7 (x,Ci))dC2^(Ci)dCi 
r (y - 7(x, Ci); s + 7(^1 Ci)) MCiXi, 

where the second equality above again follows from Lemma [2.34 - (2) . Finally, the representations in (|35[) can be 
deduced for e small enough from the relationships ([69 |) - ([70| . 
(2) Boundedness of d y lZ l (w;x,y): 

Analyzing the remainder terms lZ 2 (x;y), TZ^(x;y), lZ\(x;y), and TZ® (w; x, y), it transpires that it suffices to show 
that d y L 2 Ho(w; x + y), d y L e Ho(w;x + y), d y L £ Hi(w; x + y), and d y L e H 2 (w;x + y) exist and are uniformly bounded 
in w and y. From the definition of L £ in (|20p. one can see that, for any function H(w;y) : M 2 — > R in C£°(M. 2 ), 
d y (L e H(w;y)) exists and 



d y (L £ H(w;y)) = L e (d y H)(w;y), sup \d y L e H(w; y)\ < 00. 



w,y 

As in the proof of Lemma l2 . 31 - (4) and the relationship (|()2"j) . one can verify that Hq(w; x+y), Hi(w; x+y), H2{w; x+y) 
are C£° functions. 

In order to show that d y lZ\(w;x,y) and d y TZ^(w;x,y) exist and are bounded, it suffices that the remainder term 
TZt{z; •&) of ([39f is differentiable with respect to d and d^TZt(z; $) is bounded. The remainder term is defined as in 
(|47p , which in turn is defined as the limit as 5 — > of each of the four terms in (|4"5)) . We will show that the limit of 
the second term is indeed differentiable with respect to $ and its derivative is bounded. The other three terms can 
be dealt with similarly. As shown in (|52|) . the limit of second term in (|48 | turns out to be J f ' (i?;z,e)+/ t ' +&;z,e) 
and, thus, it suffices to show that the functions I\ ' (w;z,e) and I 4 ' (w;z,e) defined in (|49|) arc differentiable 
with respect to w and that their respective derivatives are bounded. The latter two conditions will follow from 
Lemma TB.2I together with the same arguments leading to (|50p - (|5"Tj) . ■ 
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D Proofs of other lemmas and additional needed results 



Proof of Lemma 12.31 In light of (C2-b) and the continuity of (2, £) -> £^7(2, £), we have that either ^7(2, £) < 
—5 or £^7(2, > <$, for all (2,£). Without loss of generality, we assume throughout the proof the latter case (the 
other case can be treated in the same way). Since £^7(2, £) > the mappings ( — » 7(2, Q and £ — > l( z >0 are 
both strictly increasing for all 2 € R and their inverse exists. By the same argument, we can assume without loss 
of generality that 1 + d x j(x,Q > S for all x, £ and, thus, the mapping 2 — > z + 7(2, C) =: u is strictly increasing 
and its inverse 7(w, C) exists. We now show the different assertions of the lemma: 

(1) Let us just consider 7 _1 (z, C) (the verification for 7 _1 (z, C) g° es along the same lines). Since 7(2, £) = C(z, () + z 
is continuously diffcrcntiable and condition (C2-b) is satisfied, the Inverse Function Theorem implies 7 _1 (2, C) is 
diffcrcntiablc in £ and 



^(*>o = (|| (^r^c))" = (^('.r^o))" • (76) 



Moreover, since 7 is C°°(R x R), the Implicit Function Theorem implies that 7 1 (z, £) is C°°(R x R). By implicitly 
differentiating the equation 7 _1 (z, 7(2, C)) = C with respect to 2 and then using (|76p . we get 

^7 -1 ^ a\ ^7 / --1/ ( ( ~-i 



(*.o - kr^o) ( ^ M^M) ) (77) 



(£ (*,r l (*,o) + 1) (f^rW)) 



By formally differentiating the expressions ([76]) and (|77j). it follows that 7(2, C) is C^ 1 - Indeed, by induction, one 
can verify that for any k, £ > with fc + £ > 1, 



(*> ^ 0) = (g (*, C)) " ' * ( (z, C) : 1 < k' + £' < k + 1 



where $ ( JpTg^i 7 " ( z ' C) : k! + £' < k + £, k' + £' > lj represents a polynomial function of partial derivatives of 7 of 
order greater than but no greater than k + £. Since 7(2, C) £ C^ 1 and \d^{z, £)| = ^7(2, C)l > <^ we conclude 
that 7 -1 (z,C) e C^ 1 as well. 

(2) Since £ — > 7(2, C) is assumed to be strictly increasing, 



F(j(z,j s )>o = w(J s >r 1 (z,o) = / /i e («)d« = F e (r 1 (^0). (78) 

where F e (ii) = J*J° h £ (v)dv. Clearly, the function (2, Q -> P (7(2, J e ) > Q is C£°(R x R) (being the composition of 
two C£° functions). In particular, by (f76|) , 

^p(7(^j e ) > = ^(r 1 ^, c))^(*,o = -he(r\z,o) (s.r'foo)) 1 , 

which in turn leads to the existence of T(£; z). This density admits the representation: 

f (c; z) = fteCr'c*, 0) (§2 (*, r 1 ^; o)) \ 

which is of the form (fT5|) since (^7(2, Q = d^(z, C) was assumed to be positive. Being h e 6 C£° and dj/d( G C^ 3 
and greater than J, it is clear that T e C£°. The density for 7(2, J e ) can be treated in the same was since 
£ — > 7(2, J e ) is also strictly increasing. 
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(3) Differentiating both sides of ([78)) with respect to z and using (|77j) . 

97 



9( 



Being e C£° and dj/d( £ C£° and greater than £, it is clear that ¥(j(z, J £ ) > () £ C£° . In general, we will 



have that 



(79) 



One can similarly show that (z, £) — > P (j(z, J e ) > () is C£°. 

(4) The proof goes along the same lines as that of (1) above. Note that in this case, 



(?7 1 c?7 

^ (U ' C)= 1 + (^7)(7(«,C),C)' 9C (W ' C) 



(9 C 7)(7(",0,0 

i + (a X 7)(7(«,0»0" 



(80) 



■ Proof of Lemma 13.31 Recall that the decomposition (|2"0"|) where for simplicity hereafter we shall drop the 
subindex in the operators I e and T> E . In order to ease the notation, we write Y t := X t (e, 0,x) below. Applying Ito's 
formula, we have 

f(Y t ) = f(x) + f LJ(Y u )du + f f'(Y u )a(Y u )dW u + [ I B (Y u - , C) M'(du, d(), 



where M' is the compensated Poisson measure of Z'(e) and B(y,C,) '■= fiy + liy, 0) — /(?/)■ Since f E 
and a is bounded, the second integral above is a martingale. The last integral is also a martingale. Indeed, 
io i|CI>i ^ O^u- 1 M'(du, d() is a martingale because / is bounded, while J Q J^^i & (Y u - , () M'(du, d() is a mar- 
tingale because 



B(Y U -X) h e (QdCdu< 



J\C\<1 



J\Q<U0 



due to condition (C2-a). Then, we conclude that 



\f(Y u + 7 (F„, 0P)\ 2 dp hiY u ,0\ 2 h e (0d(du < 00, 



L e fiY u )du. 



(81) 



Next, we show that L £ fiy) is bounded. Obviously, T>fiy) is bounded whenever / S and (C3) hold. Also, 



\mv)\ 



(/(y + 7(»,C))-/(i/))MCK 



< 



\f'(v + i(y,0P)W\-y(y,0\hs(0d(, 



which is clearly bounded in y due to condition (C2-a). We can then interchange expectation and integration in 
(f8Tj) . (f2"Tj) follows from the change of variables u = at. 

We now proceed to show (|2"3"|) . To this end, wc check that L £ f belongs to C\ if / G C^. Clearly, X>/ S in light 
of (C3). To show that If £C%, note that 

\d y [f(y + 7(2/, 0) - f(y)} I = + i(y, C))(i + 0„7(v, 0) - 

<ll/"lloo|7(2/,C)l + ll/'llco|a y7 (y,C)l, 



which can be bounded by C e ( || 



|/'||oo)|Cl for all ICI < e in light of (C2-a). Given that / |CMCK < 



we can interchange differentiation and integration to get 

(J/)'(y)= / (/ , (tf + 7(y,0)(i + »»7(tf,0)-/'(y))S e (0dC 



/"(y + 7(y, QP)dfa(y, + / /'(z/ + 7(j/,0H7(j/,C)MCK 
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In particular, su Py |(X/)'(t/)| < ff e (||/"||oo + Wf'W^) J |CMCK =: K and 

Sup \(L £ f)'(y)\ < ||&'||oo||/'||oo + ||&||oo||/"||oo + ||«'|U|/"||oo + |MU|/ (3) ||oo + K. 

y 

We can similarly show that sup y \(L £ f)"(y)\ < K when / £ C£, for a constant K < oo depending only on ||/( 
Halloo, ||& (i) ||oo,fori = l,...,4, and f \(\h e (QdC We also have 

L 2 J(y) = b(y)(L e fY(y) + ^l(L £ f)"(y) + J [LJ(y + 7 (y, 0) - LJ(y)] h e ({)d( 

= b(y)(L s f)'(y) + ^±( Le f)"(y) + [ [ (L £ f)'(y + 7 (y, QWMv, OMCR- 



Therefore, 

||^/||oo < \\b\\oo\\(Lefy\\oo + \\v\U\(L £ f)^oo + C4(L e fy\\ oc I \(\h e ({)d( < oo. (82) 



Since L £ f £ C 2 , by interchanging expectation and integration in (|81[) and applying again Dynkin's formula, we get 

E/(Y t ) = f(x) + J EL £ f(Y u )du = f(x) + J (L £ f(x) + J EL 2 f(Y v )dv^j du 

= f(x)+tL e f(x)+ f f EL 2 J(Y v )dvdu. 
Jo Jo 

Applying Fubini in the last integral and changing variable to a = v/t, we finally have 

E/(r t ) = f(x) + tL £ f(x) + t 2 f (1 - a)EL 2 f(Y at )da. 

Jo 



The following result is needed in order to prove Lemma IB. 21 

Lemma D.l Assume the Conditions (C1-C4) of Section^ are enforced and let <!> t : x — >• Xt(e,0,x) be the 
diffeomorphism associated with the solution of the SDE U6\). Then, for any k > 1 and compact K C M, 



lim sup sup E 

t-¥0 qEK 



d i ^7 1 

(V) 



drf 



<oo, i = 1,2. (83) 



Proof. To simplify the notation, we write X(x) = {X t (x)} t >o for {X t (e, 0, x)} t >o and fix Y s (x) := Xu_ 8 ^_(x) for 
< s < t and Yt(x) := Xq(x) = x. Wc follow a similar approach to that in the proof of Lemma 3.1 in [18] based on 
time-reversibility (see Section VI. 4 in [30] for further information). Recall that the time-reversal process of a cadag 
process V = {V s }o< s < t is given by the cadlag process 

V? = 01 s=0 + (V (t _ s) _ - V t -)l 0< . <t + (Vo - V t -)l s =f (84) 

Our main tool is Theorem VI. 4. 22 in [30j . The following notation and definitions are useful for verifying the 
assumptions in the theorem. 

Throughout, $ Sjt (-; w) :I4R denotes the diffeomorphisms defined by $ S) t(x;w) := X^ t (x;u>) where X^ t (x;u>) 
is the unique solution of the SDE 

Xl t (x) = x + f a{Xl u {x))dW u + f b{Xl u {x))du + £ 7 (A S %_ (x), AZ' U ). (85) 

Js "' s s<u<t 

The a.s. existence of this diffeomorphisms is guaranteed from (fTT]) as stated in Remark l2.2[ As usual, T s = V J\f 
and F = (.F s ) 0<s<t , where J"" = cr{W u , Z' u ;u < s} (0 < s < t) and N arc the P-null sets of J" t °. We also define 
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the backward filtration H = ( / H s ) <s<t by % s = f] s<u <t^u V cr{X t }, where (F s )o<s<t is defined analogously to 
(J-" s ) 0<s<t with W and Z' replaced with their reversal processes W* and Z' 1 . 

We are ready to show the assertions of the Lemma. First, note that, by the uniqueness of the solution of ((55)) . 
X t (x) = § Sit (X s (xj). Thus, X s (x) = <f>-}(X t (x)) G W}- 8 and, of course, X a (x) £ F s , so that a(X s (x)) e ^A^" 8 . 
Also, by Ito's formula, the quadratic covariation of W = {W s }o< s <t with a(X) := {a(X s (x))}o<s<t is given by 

a(X),w] = f a'(X u (x))a(X u (x))du = f a' (Y t _ u (x))a(Y t _ u (x))du. (86) 
J s Jo Jo 

Finally, recalling that W = {W s }o< s <t is an (F, H)-reversible semimartingale (cf. Theorem VI. 4. 20 in [30]), the 
assumptions of Theorem VI. 4. 22 in [30] are satisfied with a(X) and W in place of H and Y, respectively. By the 
theorem, we have 



a(X u (x))dW u + a(X),W = / a(X t _ u (x))dW, 



/t 
U1 



or equivalently, by ([56"]) and the change of variable v = t — u, 

1 



o 



a(X u _(x))dW u -I a'(Y v (x))a(Y v (x))dv = I <r(Y»(x))dW{. (87) 



Omitting for simplicity the dependence of the processes on x, the first term on the left-hand-side of ([87]) can be 
written as 



X.-x- f b(X u _)du- l{X u _,AZ' u ) = X (t _ s) _-X t _ + j b(X u )du+ 7(^u-,A<) 

*' 0<u<- s , ' t ~ s t-s<u<t 

= Y S -Y Q + f b(Y v )dv+ j(X (t _ v) -,AZ' t _ v ), 

where the last equality above is from the change of variable v = t — u. Then, (|87[) implies that 
Y s (x) = Y (x) - [ b(Y v (x))dv+ f a'(Y v (x))a(Y v (x))dv+ [ a(Y v (x))dWl 



l(X {t . v) -(x),AZ' t _ v ), Y Q (x)=X t -(x). 



0<v<s 

Let us write the jump component of Y in a more convenient way. To this end, note that, since Xn_ v ^-(x) + 

7(X( i _ 1 ,)- (x), AZ' t _ v ) = X t - V (x), one can express X( t _ v )-(x) in terms of the inverse j(u, () of the mapping 
z — > u := z + 7(2;, £) as follows 

Y v (x) = X (t _ v) - (x) = j(X t . v (x), AZ' t _ v ) = y(Y v - (x),AZ' t _ v ). 

Then, 

AY v (x) = 7(y fl - (x),AZ' t _ v ) - Y v - (x) = y(Y v - (x), -AZ' V ) - Y v - (x) = 7o (F„- (x), AZ' V ), 
where 7o(u, C) := l{ u ^ — C) — u an d Z' v := Z' v is the time-reversal process of {^}o<u<t- We conclude that 



Y s {x) = X t -{x) - b(Y v (x))dv+ / a'{Y v {x))a(Y v (x))dv+ / a{Y v {x))dW* v + J2 lo(Y v - (x) , AZ' V ) 



-t 

■+ 

0<v<s 



Now, define the diffeomorphism ^ s : R — > R as ^ s (il) := Y s (r)), where {V s (?7)}o<s<t is the solution of the SDE 
Y s (v)=V- f b(Y v ( V ))dv+ [ a'(Y v (, 1 ))a{Y v {i 1 ))dv+ f a(Y v (n))dW* v + £ lo(Y v -(r,), AZ' V ). 



0<v<s 
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Since, P-a.s., 

*t($t(s)) = MMx)) = *t(*t-(*)) = Yt(x) = x, for all x, 

it follows that, P-a.s., ^t(v) = *t _1 ( 7 ?) f° r au »7 G R. Furthermore, {^t(??)}t>o solves an SDE of the form (6-2) in 
[7] with their coefficients satisfying the assumptions of Lemma 10-29 therein. Finally, by Lemma 10-2-c in [7], with 
n = 2 and 5 = 1, 



sup sup E 

0<s<t r)£K 



drf 



sup sup E 

0<s<t 7]GK 



dr\ l 



= sup sup E 

0<s<t r]GK 



d l Y s {n) 



drf 



< oo, (z = 1,2). 



■ Proof of Lemma IB.21 For simplicity, we write r(0 = T(C, z) and only show the case k = 1 (the other cases 
can similarly be proved). Using the same ideas as in the proof of Proposition 1.2 in |19j . one can show that 



where 



r(C)p t (77;e,0,OdC = E(ff t fo)), 



F t (r?) : =f^^))^M»7) 



Denoting Jt(f]) :== 1 (r])/dr), note that 

^C?) = r'C^" 1 ^))^^) 2 +f(^ 1 (r 1 ))Jl(r 1 ), 
and, using (|83| and that V G C£°, it follows that sup^g^ E |i/ t '(?7)| 2 < oo. In particular, 



lim E 



H t (v + h) - H t (r)) 



= E I lim 



H t ( V + h) - H t ( V ) 



(88) 



since the set of random variables {[£?t(?7 + /i) — H t (rj)]/h : \h\ < 1} is uniformly intcgrable. Indeed, 



sup E 

\h\<l 



H t (r) + h) - H t (v) 



sup E / H' t (r) + h(3)df3 < sup E (H^rj + hf3)j < 00, 

|h|<l \Jo / W<1 

mo a] 



again due to 



can be written as 
T(0pt (m e, 0, C) dC = E (f '(Sr'fa)) (4 fa))*) + E (F^ 1 fo))3(»7) 



Then, 

d 
dr) 

It is now clear that (|41l) will hold true in light of ([83)) . ■ 

Lemma D.2 Assume the conditions (C1)-(C4) of Section^ are satisfied and letT> E andl e be the operators defined 
in i20\) . Define their respective dual operators T> E and X e as 

V 6 g(y) ■= v(y)g"(y) + (2v'(y) - b(y)) g'(y) + (v"(y) - b'(y)) g(y), 



?eg(y) 



($(7(y, 0)dvi(y,0-g(y))h e {QdC 



MO 



(g(i(y, 0) - 9(y) ff (v)(<Vr)(7(y, 0,0) YTWmMX) 



d(, 



where hereafter j(u, denotes the inverse of the mapping y — >• u := y + 7(2/, f or a fixed £ and whose existence is 
guaranteed from condition (C4). Then, the following assertions hold: 

1. D e g is well-defined and uniformly bounded for any g G and, furthermore, for any f G with compact 



support, 



g(y)V e f(y)dy = / f(y)T> E g(y)dy. 



(89) 
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2. I £ g is well-defined and uniformly bounded for any g G C'^ and, additionally, if g is integrable, then, for any 
I G CI 

' g(y)l £ f(y)dy 

3. For any g £ and f £ C\ with compact support, 

g(y)lj(y)dy= [ f'(y) [ f g( V )dr]h e (C)d(dy. (91) 

J J Jj(y,Q 

Proof. The dual relationships essentially follow from a combination of integration by parts and change of variables. 
Let us show (|9Tj) . First, note that 

g(y)l £ f(y)dy = f g(y) f f f'(y + j(y, C/?))(9 c7 )(y, WKhiO^dy 



o 



g(y)f(y + i(y, tfWa){v> ffldfciCKdCdy. 

Changing variable from y to w := 7(2/, £/3) = y + 7(1/, C/?) and applying Fubini, we get 



Finally, the second equality in (|9T|) follows from the identity 

9 € 



9 (v)dv = -9(7(w, C))da(w, C) = 0(7(10, c)) ( 9 a)(7K CU) 

1 + w)(7(w,C),C) 



Proof of Lemma lC.ll By conditioning on the times of the jumps, which arc necessarily distributed as the order 
statistics of n independent uniform [0,i] random variables, we have 

F(X t (x) >x + y\N e t =n) = — F(X t (e, {s u . . . , s„}, x) > x + y)ds n . . . ds u 

where A := {(si, . . . , s n ) : < Si < S2 < • • • < s n < t}. Hence, conditioning on T B - , 

V(X t (e, {si, s n }, x) > x + y) = E [p (x t (e, { Sl , s n }, x) > x + y\F s -) 

= E [G t - Sn {X Sn (e, {si, . . . , s„_i}, x); x, y)] , 

where Gt{z; x, y) = P (X t (e, 0, z + 7(2;, J)) > x + y). In terms of the respective densities pt(-; e, 0, C) and z) of 
X t {e, 0, C) and z + 7(2, J), we obviously have that 



G t (z;x,y) = p t {ir,e,0,Qd'nT({]z)dC= / Pt(m e, 0, C)L(C; z)d(dr). 

J J x-\-y J x-\-y J 

From Lemma [B.2| we know that there exists e small enough such that, for any S > 0, there exists B := B(e,S) < 00 
and to := to(e, 5) > such that, for all < t < to, 



sup sup I p t {ri;e,% 1 C)T{Q;z)dC < B. 

z£Mri£[x+y-8,x+y+8] J 
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The previous uniform bound for any 77 in a neighborhood of x + y and all z £ R allows us to interchange the 
differentiation and the other relevant operations (integration, expectation, etc.) so that 

d y ¥(X t (x) >x + y\N?=n) = — / d y ¥{X t {e, {a u s n }, x) > x + y)ds n ...d Sl 

1 J A 

E [d y G t _ Sn (X Sn (e, {si, . . . , s„_!}, x); x, y)] ds n ...dsx 



fn 

L J A 

nl 



and also, for any < t < t , 

Using this bound, we finally have that 



Pt-sn ( x + v'i £ > 0) C)r(C; X Sn (e, {si, . . . , s„_i}, x))d( 
\d v F(X t (x) >x + y\N£ = n)\<B. 



ds n . . . dsi 



\d y TZ t (x, y)\ < e-^* £ \d v P(X t (x) >x + y\N% = n)| { ^f- < Be^ f] < BX 3 £ t 3 . 



n— 3 



n— 3 



The proof is then complete. ■ 

Proof of Lemma 16. 11 By conditioning on the times of the jumps, which are necessarily distributed as the order 
statistics of n independent uniform [0,t] random variables, we have 



V(\X t -x\ > \ogy\N* = n) = — / P(|X t (e, {s ly . . . , s n }, x) - x\ > log y)ds n . . . ds u 

t J A 

where A := {(si, . . . , s n ) : < Si < S2 < ■ • • < s n < t}. Hence, we only need to bound 



1 



sup 



neN,te[o,i] nl Jo 
uniformly. By conditioning again, 



WX t (e,{s 1 ,...,s n },x) -x\> \ogy)dy 



X t (e, {si, . . . , s n },x) -x\> logy) =E P ( \X t (e, {s 1; . . . , s„}, x) -x\> logylj". 



< E |P [\X t . s n (e, 0, z) - x\ + | 7 (z, J)\ > logy] \ z=XMu ..., Sn _ l}>a) 
Recall the Condition (C5), we have for some constant M > and all A < 3 

supEe A|7(a: ' J)l sup < C [ e^ {x ' z) \h{z)dz < M < 00. 

X X J 

Now fix any positive constant A and t < 1, we have 

Ee |X t (e,{ Sl> ..., s „},*)-*l = P{\X t (e,{ Sl ,...,s n },x)-x\>logy}dy 

*{\X t (e,{s l7 ...,s n } 7 x) -x\> \ogy}dy 
+ / F{\X t (e,{si,...,s n },x)-x\>logy}dy 

J A 



< A 



E 



{\X t _ Sn (eAz)-x\ + | 7 (z, J)\ > logy} U^ (ei{aii ... i( , B _ l}>x) 
1 1 



dy 



< A + 2 e 5 A ? fc(1+cxp(Aie »^-- (E e x i\ x °n(e,Uu-,sn-i},x)-x\^j ^Ee Al|7(;E " 7)l ^ 

< A + 2il/e^ A i fc(1+oxp(Aie)) — - fEe Al|Xs " (e ' { 



{3i,...,s n _i},a;)— a; | 
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Above, we used (|17j) for the last inequality with A = Ai = 1 + a, where < a < 2 is to be chosen later. Now we 
iterate the above procedure by taking Ai = (1 + a) 1 , i = 1, 2, n, at each step, and choose A„ = (1 + a) n = e. We 
conclude that there exists a large enough constant C independent of n and t such that 

/•oo 

/ ¥{\X t (£,{si,...,s n },x)-x\ > \ogy}dy < C 
Jo 

In what follows, we only need to show C n (l/a) n /n! — ► as n —> oo. Recall that a = e 1 /™ — 1. We have 

~ n ( C + log — ] as n — > oo . 
\ n J 

On the other hand, we know log??.! ~ n 2 /2 as n — > oo. The proof is then complete. ■ 

References 

[1] Y. Ait-Sahalia and J. Jacod. Testing for jumps in a discretely observed process. Annals of Statistics, 37:184-222, 
2009. 

[2] Y. Ait-Sahalia and J. Jacod. Is Brownian motion necessary to model high frequency data? Annals of Statistics, 
38:3093-3128, 2010. 

[3] D. Applcbaum, Levy Processes and Stochastic Calculus, Cambridge University Press, 2004. 

[4] O.E. Barndorff-Nielsen and N. Shcphard. Econometrics of testing for jumps in financial economics using 
bipower variation. Journal of Financial Econometrics, 4(l):l-30, 2006. 

[5] H. Berestycki, J. Busca and I. Florent, Asymptotics and calibration of local volatility models, Quantitative 
Finance, 2:61-69, 2002. 

[6] H. Berestyki, J. Busca and I. Florent, Computing the Implied Volatility in Stochastic Volatility models, 
Communications on Pure and Applied Mathematics, Vol. LVIL1352-1373, 2004. 

[7] K. Bichteler, J.B. Gravereaux and J. Jacod: Malliavin Calculus for Processes with Jumps, Stochastics Mono- 
graphs Volume 2. Gordon and Breach Science Publishers (1987). 

[8] P. Carr and L. Wu, What type of process underlies options? A simple robust test, Journal of Finance, 
58(6):2581-2610, 2003. 

[9] R. Cont and P. Tankov, Financial Modelling with Jump Processes, Chapman & Hall, 2004. 

[10] J. Feng, M. Forde, and J. P. Fouque, Short maturity asymptotics for a fast mean reverting Heston stochastic 
volatility model, SIAM Journal on Financial Mathematics, 1:126-141, 2010. 

[11] J.E. Figueroa-Lopez and M. Forde, The small- maturity smile for exponential Levy models, Submitted, Available 
at arXiv:1105.3180vl [q-fin], 2011. 

[12] J. Figueroa-Lopez, and C. Houdre, Small-time expansions for the transition distribution of Levy processes, 
Stochastic Processes and their Applications, 119:3862-3889, 2009. 

[13] J.E. Figueroa-Lopez, R. Gong and C. Houdre, Small-time expansion of the distributions, densities, and option 
prices of stochastic volatility models with Levy jumps, Stochastic Processes and their Applications, 122:1808- 
1839, 2012. 

[14] Forde, M., and Jacquicr, Small-time asymptotics for an uncorrelated local-stochastic volatility model, Appl. 
Math. Finance, 18, 517-535, 2011. 




los 



34 



[15] Gao, K., and R. Lee, Asymptotics of Implied Volatility in Extreme Regimes, Preprint, Available at http:// 
papers.ssrn.com/, 2011. 

[16] J. Gatheral, E. Hsu, P. Laurence, C. Ouyang and T-H. Wang, Asymptotics of implied volatility in local 
volatility models, To appear in Mathematical Finance, 2009. 

[17] Y. Ishikawa, Asymptotic behavior of the transitions density for jump type processes in small time, Tohoku 
Math. J. 46:443-456, 1994. 

[18] Y. Ishikawa, Density estimate in small time or jump processes with singular Levy measures, Tohoku Math. J. 
53:183-202, 2001. 

[19] R. Leandre, Densite en temps petit d'un processus de sant, In: Seminaire de Probabilitts XXI, Lecture Notes 
in Math. (J. Azema, P. A. Meyer, and M. Yor, eds.), Springer, Berlin, 1247:81-99, 1987. 

[20] J. P. Lepeltier and R. Marchal Problemes de martingales associees a un operateur integro diffcrcnticl, Ann. 
Inst. Henri Poincare Probab. Stat. B. 12:43-103, 1976. 

[21] S. Levendorskii, American and European options in multi-factor jump difussion model, near expiry, Finance 
Stoch., 12:541-560, 2008. 

[22] P. Marchal. Small time expansions for transition probabilities of some Levy processes. Electronic communica- 
tions in probability, 14:132-142, 2009. 

[23] A. Medvedev and O. Scailllet, Approximation and calibration of short-term implied volatility under jump- 
diffusion stochastic volatility, The review of Financial Studies, 20(2):427-459, 2007. 

[24] J. Muhle-Karbe and M. Nutz, Small-time asymptotics of option prices and first absolute moments, Preprint, 
2010. 

[25] D. Nualart. The Malliavin Calculus and Related Topics. Springer- Verlag, New York, USA, 1995. 

[26] B. Okscndal and A. Sulem: Applied Stochastic Control of Jump Diffusions, Springer (2005). 

[27] J. Picard, On the existence of smooth densities for jump processes, Probab. Theory Related Fields, 105:481-511, 
1996. 

[28] J. Picard, Density in small time at accessible points for jump processes, Stochastic processes and their 
applications, 67:251-279, 1996. 

[29] M. Podolskij. New Theory on estimation of integrated volatility with applications. PhD thesis, Ruhr-Universitat 
Bochum, April 2006. 

[30] P. Protter. Stochastic Integration and Differential Equations. Springer, 2004. 2nd Edition. 

[31] M. Roper, Implied volatility: small time to expiry asymptotics in exponential Levy models. Thesis, University 
of New South Wales, 2009. 

[32] L. Riischendorf and J. Woerner. Expansion of transition distributions of Levy processes in small time. Bernoulli, 
8:81-96, 2002. 

[33] I. Shigekawa: Stochastic Analysis, Translations of Mathematical Monographs Volume 224. AMS (2004). 

[34] P. Tankov, Pricing and hedging in exponential Levy models: review of recent results, To appear in the 
Paris-Princeton Lecture Notes in Mathematical Finance, Springer 2010. 

[35] J. Yu, Closed-form likelihood approximation and estimation of jump-diffusions with an application to the 
realignment risk of the Chinese Yuan. Journal of Econometrics, 141:1245-1280, 2007. 



35 



